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ABSTRACT 

In this paper we present results from an IRAM Plateau de Bure millimetre-wave Interferom¬ 
eter (PdBI) survey for CO emission towards radio-detected submillimetre galaxies (SMGs) 
with known optical and near-infrared spectroscopic redshifts. Five sources in the redshift 
range z ~ 1-3.5 were detected, nearly doubling the number of SMGs detected in CO. 
We summarise the properties of all 12 CO-detected SMGs, as well as 6 sources not de¬ 
tected in CO by our survey, and use this sample to explore the bulk physical properties of 
the SMG population as a whole. The median CO line luminosity of the SMGs is (Lpo) = 
(3.8 ± 2.0) X 10^°Kkms“^pc^. Using a CO-to-H 2 conversion factor appropriate for star- 
burst galaxies, this corresponds to a molecular gas mass (M(H 2 )) = (3.0 ± 1.6) x 10^° Mq 
within a ~ 2 kpc radius, about four times greater than the most luminous local ultraluminous 
infrared galaxies (ULIRGs) but comparable to that of the most extreme high-redshift radio 
galaxies and QSOs. The median CO FWHM linewidth is broad, (fwhm) = 780±320 km s“^, 
and the SMGs often have double peaked line profiles, indicative of either a merger or a disk. 
From their median gas reservoirs (~ 3 x 10^° Mq ) and star-formation rates (> 700 Mq yr“^) 
we estimate a lower limit on the typical gas-depletion time scale of > 40Myr in SMGs. This 
is marginally below the typical age expected for the starbursts in SMGs, and suggests that 
negative feedback processes may play an important role in prolonging the gas consumption 
time scale. We find a statistically-significant correlation between the far-infrared and CO lu¬ 
minosities of the SMGs which extends the observed correlation for local ULIRGs to higher 
luminosities and higher redshifts. The non-linear nature of the correlation implies that SMGs 
have higher far-infrared to CO luminosity ratios, and possibly higher star-formation efficien¬ 
cies, than local ULIRGs. Assuming a typical CO source diameter of 6* 075 (D ^ 4kpc), 

we estimate a median dynamical mass of {Mdyn) — (1-2 ± 1.5) x 10^^ Mq for the SMG 
sample. Both the total gas and stellar masses imply that SMGs are very massive systems, 
dominated by baryons in their central regions. The baryonic and dynamical properties of these 
systems mirror those of local giant ellipticals and are consistent with numerical simulations 
of the formation of the most massive galaxies. We have been able to impose a lower limit of 
> 5 X 10“® Mpc“^ to the co-moving number density of massive galaxies in the redshift range 
z ^ 2-3.5, which is in agreement with results from recent spectroscopic surveys and the most 
recent model predictions. 

Key words: galaxies: starburst - galaxies: formation - cosmology: observations - cosmol¬ 
ogy: early Universe 
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1 INTRODUCTION 

The discovery of extragalactic carbon monoxide (CO) rota¬ 
tional line emission (Rickard et al. 1975) and, in particular, the 
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first detections of CO at cosmologically significant redshifts in 
IRASF10214+4724 at 2 = 2.29 (Brown & Vanden Bout 1991; 
Solomon, Downes & Radford 1992a), the Cloverleaf at 2 = 2.56 
(Barvainis et al. 1994) and BRl 1202—0725 at 2 = 4.69 (Omont et 
al. 1996) revealed the potential of CO as a tracer of molecular gas 
in the early Universe. Since those pioneering efforts, progress has 
been slow due to the severe observational obstacles: the faintness 
of the CO emission, except for cases were the source is gravita¬ 
tionally magnified; inaccurate spectroscopic redshifts coupled with 
the small instantaneous bandwidth of most modern-day correla¬ 
tors; and in some cases, an unfortunate combination of redshift and 
available receiver coverage, as well as atmospheric transparency, at 
millimetre wavelengths. As a result only 30-or-so 2 > 1 objects 
have been detected to date, most of which have been extremely lu¬ 
minous, often gravitationally lensed, QSOs (e.g. Omont et al. 1996; 
Guilloteau et al. 1997, 1999; Downes et al. 1999; Cox et al. 2002; 
Bertoldi et al. 2003; Beelen et al. 2004) and high-redshift radio 
galaxies (HzRGs — e.g. Papadopoulos et al. 2000; De Breuck et 
al. 2003a). 

The slow increase in the number of CO detections contrasts 
with the rapid growth in samples of high-redshift galaxies selected 
through continuum observations at submillimetre (submm) wave¬ 
lengths, which detect thermal emission from dust. The advent of 
large-format submm/mm cameras, SCUBA (Holland et al. 1999) 
and MAMBO (Kreysa et al. 1998), revealed the presence of a 
significant population of dust-enshrouded, and therefore hitherto 
undetected, galaxies at high redshifts (e.g. Smail, Ivison & Blain 
1997; Barger, Cowie & Sanders 1999; Bertoldi et al. 2000). Today 
such observations are routine, and several hundred submm-selected 
(or SCUBA) galaxies (SMGs) have been discovered (e.g. Blain et 
al. 2002; Scott et al. 2002; Webb et al. 2003; Borys et al. 2003). 

Until recently, only two SMGs had been detected in CO 
(Prayer et al. 1998, 1999; Downes & Solomon 2003; Genzel et 
al. 2003), largely owing to the extreme faintness of SMGs in the 
optical and the difficulties in obtaining reliably spectroscopic red¬ 
shifts. However, a major step forward was made by Chapman et al. 
(2003b, 2005) who used the highly-efficient, blue-sensitive LRIS- 
B spectrograph on the Keck Telescope to obtain spectroscopic red¬ 
shifts for a large sample of SMGs. This motivated a major sur¬ 
vey at the IRAM Plateau de Bure Interferometer (PdBl) to look 
for CO emission from this high-redshift sample. The observational 
cycle involves: identifying a robust radio counterpart to an SMG 
in deep VLA radio maps (see e.g. Ivison et al. 2002); placing a 
LRIS-B/Keck slit on this position to obtain a spectroscopic redshift 
(Chapman et al. 2003b, 2005); frequently confirming the redshift in 
the near-infrared, usually via redshifted Ha (Simpson et al. 2004; 
Swinbank et al. 2004); and, finally, searching for redshifted CO 
emission with PdBI. The initial results from this CO programme 
were described in Neri et al. (2003). Although an expensive pro¬ 
cess in terms of telescope time, this is currently the only feasible 
and effective route for detecting CO from SMGs. 

CO observations can potentially provide unique information 
about the enigmatic SMG population. First and foremost, CO traces 
the bulk of the molecular gas in SMGs: the high-level (J > 2) tran¬ 
sitions arise in the warm and dense gas while lower level J lines 
probe the quiescent and likely cooler gas which may lurk in the 
outskirts of SMGs. Detecting and mapping this molecular emis¬ 
sion enables us to precisely determine the spatial and kinematic 
location of the gas-rich components within an SMG. About two- 
thirds of SMGs are found to be large, morphologically complex 
systems in the optical/near-infrared and/or the radio, with typically 
one or more companions (e.g. Smail et al. 1999; Ivison et al. 2002; 


Chapman et al. 2003c, 2004). For example, in one SMG, Neri et 
al. (2003) identified CO emission coincident with a second, fainter 
radio source ^ 4” away from the radio counterpart for which the 
spectroscopic redshift had been found. Secondly, CO observations 
yield fairly accurate estimates of the amount of molecular gas avail¬ 
able to fuel the starburst and/or the AGN responsible for the large 
far-infrared luminosities. While the conversion factor between the 
CO luminosity and molecular gas mass is uncertain, CO observa¬ 
tions clearly provide a much better constraint on the amount of gas 
present in SMGs than the estimates based on submm continuum 
observations and an adopted spectral energy distribution (SED) and 
gas-to-dust ratio. 

The first two CO detections of SMGs revealed the presence 
of copious amount of molecular gas (^ 10^° M© ), suggesting that 
intense star formation is occuring in these systems (Frayer et al. 
1998, 1999). Furthermore, from a reliable estimate of the gas reser¬ 
voir in an SMG we can say something about the gas exhaustion 
time-scale, i.e. the duration of the submm-luminous phase. This, in 
turn, allows us to make an educated guess about the possible de¬ 
scendants of SMGs, and thus place them in an evolutionary context 
with other high- and low-redshift galaxy populations. 

Observations of the shape and width of CO lines also pro¬ 
vides important information about the kinematics in SMGs (Neri 
et al. 2003; Tacconi et al. 2005). In particular, if the CO emis¬ 
sion is spatially resolved we can use that, in conjunction with 
the width of the line profile, to constrain the dynamical mass of 
the host galaxy. Estimates of the dynamical mass based on CO 
are likely to be ‘cleaner’ than estimates from optical/near-infrared 
spectroscopy, which are prone to extinction by dust and the effects 
of non-gravitational motions in the emitting gas, such as outflows. 
The best example of resolved CO emission in a SMG (Genzel et 
al. 2003) shows gas extended on scales of ~ 3-5 kpc and that 
most of the dynamical mass (^ 3 x 10^^ M© ) is baryonic. Such 
large, widely-distributed gas reservoirs suggest that the brightest 
SMGs are not merely high-redshift replicas of the local population 
of ULIRGs. However, a representative picture of the gas distribu¬ 
tion in SMGs will have to await high-resolution CO observations 
of a large sample of SMGs. 

Since CO observations provide a means of ‘weighing’ galax¬ 
ies at high redshifts, both in terms of their baryonic gas mass con¬ 
tent and their total dynamical mass, they can be used to help piece 
together a picture of the mass assembly of massive galaxies in the 
early Universe. In the classical cold-dark matter (CDM) scenario 
of structure formation (White & Frenk 1991), massive spheroidal 
galaxies are the end products of a gradual build-up of mass via 
merging and, as a result, form late in the history of the Universe. 
However, as pointed out by Genzel et al. (2003), at least some 
SMGs appear to be massive, baryon-dominated galaxies already at 
2 ; ^ 2-3. The discovery of massive, evolved early-type galaxies at 
z Z 1-5 (e.g. Cimatti et al. 2004), which could be the descendants 
of z > 2 SMGs, suggests that some massive spheroids were not 
formed at a ^ 1. In fact, if the brightest 25 per cent of the SMGs 
are massive baryonic systems, then their abundance indicates that 
the build up of massive galaxies in the early Universe was much 
faster than previously expected, and could pose a serious challenge 
for current models of gas processing in galaxy formation (e.g. Cole 
et al. 1994; Kauffmann et al. 1999; Baugh et al. 2004). 

In this paper we present the most recent results from a system¬ 
atic survey of CO emission towards radio-identified SMGs with 
spectroscopic redshifts in the range 2 ; ^ 1-3.5 (corresponding to 
40-10% of the age of the Universe at the respective epochs), and 
use these to address the issues outlined above. This is part of a 
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major effort currently being undertaken at the PdBI with the aim 
of detecting and imaging CO emission towards ~ 20-30 SMGs. 
In !j2|and|3 we describe the observations and derive properties for 
each new source observed. In !|4|we define a sample of 12 sources 
consisting of all SMGs detected in CO to date, from which the av¬ 
erage physical properties of the bright SMG population are derived. 
We discuss their properties and compare with those of other galaxy 
populations in m Finally, m discusses the impact of our obser¬ 
vations on models of galaxy formation and evolution. Throughout 
we adopt a flat cosmology, with ~ 0.27, Ha = 0.73 and 
Ho = 71 kms“^ Mpc~^ (Spergel et al. 2003). 


2 OBSERVATIONS 

In total, 11 SMGs were targeted for CO observations with PdBI 
(see Table 0, drawn from five independent submm surveys: the 
SCUBA Lens Survey (SMM102396-0134 — Small, Ivison 
& Blain 1997); the Hubble Deep Field (SMM J12360-F6210 
— Chapman et al. 2003a); the Hawaii Survey Fields, SSA13 
and SSA22 (SMM J13120-F4242, SMM J13123-F4239 — 
Chapman et al. 2005; SMM J22174+0015 — Barger, Cowie 
& Sanders 1999); The SCUBA UK 8mJy Survey of the 
Lockman Hole East and ELAISN2 (SMM J10523-F5722, 
SMMJ10524-F5724, SMM J16363+4055, SMMJ16363-F4056, 
and SMMJ16366+4105 — Scott et al. 2002); and a 1200-Atm 
MAMBO survey of the same two fields (SMM 116371+4053^ — 
Greve et al. 2004b). 

Deep radio imaging was used to identify radio counterparts 
to all the targeted SMGs and accurately locate these relative to 
optical/near-infrared reference fields (Small et al. 1999; Chapman 
et al. 2003a; Ivison et al. 2002). Subsequent spectroscopy with 
Keck/LRIS-B (Oke et al. 1995) and CEHT/OSIS-V in the case 
of SMM J02396—0134 (Soucail et al. 1999) provided UV spectro¬ 
scopic redshifts for all sources (see Chapman et al. 2003b, 2005). A 
subset of the sample were spectroscopically observed in the near- 
infrared to provide more reliable systemic redshifts from the wave¬ 
lengths of redshifted Hq or [Olll] A5007A emission lines (Swin- 
bank et al. 2004) to aid in our search for redshifted CO emission. 
Here, and for the remainder of this paper, we take the systemic red- 
shift to be equivalent to the CO redshift. 

From our sample of 11 sources, 10 were part of our 
statistically-complete CO survey and were targeted during two ob¬ 
serving campaigns, during the winter of 2002-03 and the summer 
of 2003, in good to excellent weather. The observations were done 
in D configuration in order to maximise the sensitivity, and used 
five of the six available antennae, giving a total of 10 baselines. We 
stress that the survey is on-going and that the 10 sources reported 
here are merely those which were observed during the 2003 season. 
In addition, we included an 11th source, SMM J02396—0134, for 
which data was originally obtained by J.-P. Kneib and G. Soucail 
during summer/autumn of 1999, and these were retrieved from the 
PdBI data archive. Observations of this source were done in D and 
C configuration. The details of the observations of all 11 sources 
are summarised in TableQ 

While our main goal was to detect redshifted CO emission in 
the 3-mm waveband, we also used the 1.3-mm receivers to attempt 
to determine accurate continuum positions and fluxes or, in cases 
where a higher-J CO line would coincide with the 1.3-mm band, 

^ This source was denoted MM J16371+4053 in Greve et al. (2004b) 


Table 1. Log of the PdBI observations for the 11 SMGs analysed in this 
paper. The on-source integration time {tint) is the observing time for the 
equivalent six element aiTay. 


Source 

Observing Dates 

tint 

(hr) 

Detection? 

SMMJ02396-0134 

1999 Jun 20, 24-26, 29, Aug 24, 27, Sep 7 

lA 

Y 

SMMJ10523+5722 

2003 Apr 25, May 06 

16.6 

N 

SMMJ10524+5724 

2003 Mar 27 

3.9 

N 

SMMJ12360+6210 

2003 Jun 01, 03, 18, Aug 30, Sep 01 

9.4 

N 

SMMJ13120+4242 

2003 May 17, 22-23,26 

9.2 

Y 

SMMJ13123+4239 

2003 Jun 23, 26, 27 

4.1 

N 

SMMJ16363+4055 

2003 Mai-31, Apr 03,21 

11.2 

N 

SMMJ16363+4056 

2003 Aug 11,22, 24 

9.6 

N 

SMMJ16366+4105 

2003 Apr 13-14, 23, May 1 

14.9 

Y 

SMMJ16371+4053 

2003 Jul 25-26, 29 Aug 1,2, Sep 4 

17.1 

Y 

SMM J22174+0015 

2003 May 9-13, 29-30, Jun 4, Sep 12, 14 

17.3 

Y 


to search foremission from this transition, e.g. SMM J04431+0210 
(Neri et al. 2003). To achieve our goals the correlators were config¬ 
ured for line and continuum observations and simultaneously cov¬ 
ered 580 MHz in the 3-mm and 1.3-mm bands: this corresponds to a 
typical velocity coverage of 1700 and 750 km s“^ at 3 and 1.3 mm, 
respectively. 

Where a near-infrared spectroscopic redshift was unavailable 
we have to correct for the likely systematic blueshifts of UV line 
features relative to the CO emission. We therefore tuned the 3-mm 
receivers slightly redward of the measured spectroscopic redshift. 
This meant that for a typical redshift uncertainty of A 2 = 0.005, 
any source at 2 > 1 would have its CO line peak fall within 0.5 x 
580 MHz = 290 MHz of the 3-mm band centre. In a few instances, 
a line was detected at the edge of the band pass. The frequency 
setting was then adjusted to centre the line in the bandpass, and 
the source was re-observed. A source was typically observed for 
2-3 tracks (10-18 hrs). If no signal had been detected after this, the 
source was not pursued further. 

All data reduction employed the IRAM GILDAS software 
(Guilloteau & Lucas 2000). This involved careful monitoring of the 
quality of the data throughout a track, and subsequent flagging of 
any bad and high phase-noise visibilities. For passband calibration 
we typically used one or more bright quasars. Phase and amplitude 
variations within each track were calibrated out by inter-leaving 
reference observations of nearby fainter quasars every 20 m. In the 
best conditions, the typical rms phase noise per baseline at 3 mm 
was < 10°, increasing to < 40° in the worst cases. Observations of 
the primary calibrators, 3C454.3, 3C345, 3C273, and MWC 349, 
were used to determine the flux scale. Finally, naturally weighted 
data cubes were created using GILDAS. 


3 THE SOURCES 

No sources were detected at L3-mm, and the data were of so poor 
quality that even upper limits were useless. For the remainder of 
this paper we shall therefore only discuss the 3-mm data. For each 
source the continuum level at 3-mm was estimated by extrapolating 
from its submm flux and assuming a spectral index of /3 = 1.5. In 
all cases we find the contribution from the continuum emission to 
the CO line at 3-mm to be less than < 5 per cent, and therefore 
negligible. Of the 11 sources whose observations are analysed in 
this paper, CO emission was reliably detected in five. Including the 
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first three (out of three) detections from our survey from Neri et al. 
(2003) yields a detection rate of 57% (8/14) to date. 

The failure to detect CO in some of our sources is unlikely 
to be due to the optical redshift being wrong, since the quality of 
the spectra is generally high (Chapman et al. 2005). Similarly, the 
possibility that the identification is wrong, and therefore the red- 
shift is for the wrong object, appears to be remote given the excel¬ 
lent correspondence between the radio and submm (e.g. Ivison et 
al. 2002), and the extremely small likelihood of finding a 2 : ~ 2 
galaxy at that position, which is not related to the submm emis¬ 
sion. Instead we feel the explanation for the missing CO emission 
is due to either too large a velocity offset between the CO and op¬ 
tical emission (putting the former outside our selected correlator 
coverage, see 114.2> . a CO line with a width comparable to our cor¬ 
relator bandwidth, or simply that the sources are too faint in CO to 
be detected in the given integration time. In the latter case, the 6 
non-detections can provide useful upper limits on their CO lumi¬ 
nosity and gas masses and as a result we have included them in the 
analysis in this paper. 

In remainder of this section we shall briefly describe the prop¬ 
erties of the five SMGs for which new CO detections are presented 
in this paper as well as the 6 non-detections. 

SMMJ02396-0134 The CO(2-l) spectrum of 
SMM J02396—0134 is shown in Fig. 0. The blue wing of 
the line is missing due to the limited correlator bandwidth available 
in 1999, but nonetheless it is evident that the line is double peaked. 
The CO redshift, zco = 1.062 ± 0.002, was defined as the flux- 
weighted redshift, see section l4~^ for details. The widths of the two 
peaks were estimated by simultaneously fitting two Gaussians to 
the spectrum, which yielded FWHM values 180 ± 20 kms“^ and 
430 ± 50 km s“ ^ for the blue-shifted and red-shifted peaks, respec¬ 
tively. The velocity offsets of the two peaks were —280 km s“^ and 
90 kms“^ , respectively, measured relative to the CO redshift. 
From the double Gaussian fit the velocity-integrated CO(2— 1) line 
flux is Ico = 3.4 ± 0.3 Jy kms“^. A single Gaussian fit to the line 
profile leads to a FWHM value of 780 ± 60 kms“^ . Due to the 
truncated line we can only impose a lower limit ~ 870 kms“^ on 
the full width at zero intensity (FWZI). 

The emission integrated over the line from —460 kms“^ to 
110kms“^ is shown as contours overlaid on a 7T-band im¬ 
age in Fig. HJi. The emission is detected at >10cr, where cr = 
0.40mJybeam“^, and coincides with the optical/near-infrared 
counterpart to within the relative astrometrical errors. The source 
appears extended in the north-south direction, and an elliptical 
Gaussian fit to the CO emission in the image plane yields a source 
size of 9'/5 x 5'/5, which is slightly larger than the 8'.'3 x 5'.'3 syn¬ 
thesized beam. The two velocity components of the source have 
a slight spatial offset, with the blue-shifted peak (integrated from 
—460 to —120 km s“^ ) 1.2" north of the centroid of the integrated 
line emission and the red-shifted peak (integrated from —120 to 
110 km s“^ ) offset 0.3" east and 1.6" south of the central position. 
The positional uncertainties in the maps are 0.3-0.4", suggesting 
that the observed offsets are significant. 

SMMJ13120+4242 The CO(4-3) spectrum of 
SMM J131204-4242 is shown in Fig. m- A single Gaussian pro¬ 
vides an excellent fit to the data, as measured by the statistic, and 
yields a velocity-integrated flux density, 1.7 4z 0.3 Jy kms“^ and a 
line width of FWHM ~ 530 ± 50 kms“^ . We find a systemic CO 
redshift of Zco = 3.408 4= 0.002. 

The velocity-integrated CO(4—3) emission yields a > 6-a de¬ 
tection (a — 0.32 mjy beam“^) at position which is offset 0.8" to 
the north-east of the optical counterpart (Fig.|2>)- The positional 


error is ~ 1", and so the offset is not significant. The CO emission 
is unresolved. 

SMMJ163664-4105 The CO(3-2) spectrum of 
SMM 1163664-4105 (Fig.0) is very broad with the blue edge of 
the line ending abruptly. The spectrum exhibits a double-peaked 
line similar to SMM 02396—0134, although less distinct. The 
CO redshift is z = 2.450 4z 0.002. A simultaneous fit of two 
Gaussians provides a superior fit to the spectrum than a single 
Gaussian, and places the blue-shifted and red-shifted peaks at 
velocity offsets of —270kms“^ and 4-210kms“^, respectively. 
The velocity-integrated line flux is Ico = 1.8 ± 0.3 Jy kms“^ . 
A single Gaussian fitted to the spectrum yields a line width 
FWHM=870 4r 80 kms“^ - the second largest line width of the 
sample. 

The integrated CO emission, shown in Fig. HI;, is detected at 
the > 8-a level {a = 0.21 mjybeam“^), and coincides with the 
radio and faint optical/near-infrared counterpart. We imaged the 
blue-shifted and red-shifted peaks separately, using different veloc¬ 
ity cut-offs, in order to look for positional offsets between the two, 
but found no evidence for a velocity gradient across the source. It 
is possible that part of the CO emission is related to the faint near- 
infrared source ~ 2" south-west of the CO centroid, the compact 
radio counterpart displays a similar faint tail of emission in this 
direction. SMM J16366+4105 is a prime target for further, higher- 
resolution CO observations. 

SMMJ16371+4053 The detection of CO(3-2) towards 
SMM J163714-4053 is the first for a MAMBO-selected source. 
The CO(3—2) line profile is shown in Fig.Q. The flux-weighted, 
systemic CO redshift is Zco ~ 2.380 ± 0.004. While clearly de¬ 
tected, the poor signal-to-noise ratio prevents us from defining 
the line profile, except to say that the line appears broad with 
FWZI> 900kms“^. A Gaussian fit yields a formal line width 
of FWHM~ 830kms“^ and a velocity-integrated line flux Ico = 
1.0±0.2Jy kms"i . 

The integrated emission is shown in Fig. m The source is 
detected at > 5cr (ct = 0.18mJybeam“^) and coincides with the 
optical/near-infrared counterpart to within the positional errors. 

SMM J22174+0015 The CO(3—2) spectrum of this source 
shown in Fig. 0 is the weakest of our detections. After heavily 
binning the spectrum, we estimate the line to be significant at the 
> 4 -(t level. We measure the flux-weighted, systemic CO redshift 
as Zco = 3.099 ± 0.004. A Gaussian fit to the spectrum yields a 
line width of FWHM = 780 db 100 kms“^ and a line flux density 
of Ico = 0.8 ± 0.2 Jy kms“^ . 

Contours of the velocity-integrated CO emission are shown 
on a TT-band image in Fig. Ht. The CO emission is detected at 
the ~ 4.7-0" level (a = 0.17mjy beam“^), at a position offset by 
~ 2 . 5 " north-west from the optical counterpart, but only ~ 0.8" 
away from the radio position (phase centre). 

Non-detections In Table 2 we list the 3o upper limits on 
the 3-mm CO line flux spatially coincident with the radio posi¬ 
tion as well as the redshift range searched for each of the 6 non¬ 
detections. The upper limits were calculated using the equation 
Ico = 3a (JuAufwhm)^^^, where a is the channel-to-channel rms 
noise, Sv the velocity resolution and Aufwhm the line width (see e.g. 
Seaquist, Ivison & Hall 1995). In our calculation of the upper limits 
we used the rms noise in the spectra binned to 100 MHz resolution 
and adopted a line width of 500 km s“^. The integrated CO (3—2) 
map of SMM J12360+6210 shows a ~ 4 (t peak slightly eastwards 
of the phase-centre, and is thus possibly another tentative detection, 
albeit marginal. In the remaining cases no line emission is observed 
at a significance > 3cr. 
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Figure 1. CO line profiles of the 12 SMGs detected in CO to date. The profiles have been plotted on similar flux- and velocity-scales in order to facilitate an 
easy comparison between individual sources, to achieve this the spectra of SMMJ14011-1-0252 and SMM J22174-I-0015 have been scaled by a factor xl/3 
and X 3, respectively. The LSR velocity-scale is relative to the CO redshift. Optical redshifts are denoted by arrows and the horizontal bars on top of the arrows 
indicate the uncertainty of the redshifts. The redshifts derived from the CO lines are given in Table 2. The solid curves represent the best-fits to the spectra 
using either single or double Gaussian profiles. The spectra shown in c), d) and g) are taken from Neri et al. (2003); h) is taken from Kneib et al. (2004), while 
b), f) and k) are from Genzel et al. (2003), Downes & Solomon (2003) and Andreani et al. (2000), respectively. The spectra have been binned into 20 MHz 
bins, except for f), h) and 1) where the frequency binning is 7, 10 and 40 MHz, respectively. 

© 0000 RAS, MNRAS 000, 000-000 












6 Greve et al. 


C0(2-1) SMM J02396-0134 



C0(3-2) SMM J16371+4053 



C0(4-3) SMM J13120144242 



CO(3-2) SMM J2217440015 



C0(3-2) SMM J1636644105 



Figure 2. Contours of the velocity integrated CO emission overlaid on 30” X 30” iC-band images of the five new CO-detected SMGs presented in this paper. 
The iC-band images are centred on the radio positions, while the open crosses mark the position of the optical/near-infrared counterparts for which the spec¬ 
troscopic redshifts were obtained. In all panels the contours start at 3 (t and increase in steps of Itr, where cr = 0.40, 0.32, 0.21, 0.18, and 0.17mJy beam“^ 
in panels a), b), c), d), and e), respectively. The synthesized beams are (panels a) to e), shown as hatched ellipses) 873 X 573 at position angle 15° (east of 
north), 679 X 478 at 66°, 774 X 379 at 62°, 774 X 570 at 61°, and 970 X 676 at 65°. The (Icr) uncertainties on measured positions in the optical and mm 
frames are 07l and 073, respectively. The 7f-band data for a), b) and e) were published in Smail et al. (2002), Small et al. (2004) and Chapman et al. 
(2003a), respectively, while c) and d) were published in Ivison et al. (2002). 


4 SAMPLE PROPERTIES 
4.1 The sample 

The main goal of this paper is to assemble a large sample of SMGs 
observed in CO. In order to achieve this, we add to the 11 new 
sources in Table Q from our survey (detections as well as non¬ 
detections) and the three earlier CO detections from this survey 
described in Neri et al. (2003), the three SMGs detected in CO 
prior to our survey: SMM J02399—0136 (Prayer et al. 1998; Gen- 
zel et al. 2003), SMM J1401H-0252 (Prayer et al. 1999; Downes & 
Solomon 2003) and SMM J16359+6612 (Sheth et al. 2004; Kneib 
et al. 2004). The latter is a triply-imaged lensed galaxy with CO 
detected towards all three components, and in Pig. Eh we have dis¬ 
played the combined (and therefore highest signal-to-noise) spec¬ 
trum of all three components as given in Kneib et al. (2004). Por 
the velocity-integrated line flux we have adopted the value of the 
component with the strongest detection, i.e. component B (Kneib 
et al. 2004) - see Table 2. 

We also include the extremely red object, PIR 10, at 2 : = 1.443 
(EROJ164504-4626; Hu & Ridgway 1994). Initially studied be¬ 
cause of its extremely red optical-near-infrared colours, subsequent 
observations with SCUBA and Hubble Space Telescope revealed 
that HR 10 is a powerful SMC with a distorted rest-frame UV mor¬ 
phology (e.g. Cimatti et al. 1998; Dey et al. 1999). Hence, HR 10 


bears all the characteristics of a SMC and would have been detected 
in any of the SCUBA surveys to date. The CO(5—4), CO(2—1) 
and CO(l—0) lines have all been detected towards this source (An- 
dreani et al. 2000; Greve et al. 2003). 

Pinally, we note that Hainline et al. (2004) detected CO(3—2) 
towards the submm source SMM J04135-I-10277, the first type-1 
QSO to be selected at submm wavelengths (Knudsen, van der Werf 
& Jaffe 2003), as part of an Owens Valley Millimeter Array survey 
of CO emission towards high-redshift QSOs. SMM J04135+10277 
is not only one of the brightest submm sources known but also an 
extremely luminous CO source. While SMM J04135-I-10277 was 
discovered in a blank-field (albeit lensing-assisted) submm survey 
and thus meets the selection criteria for our sample, only about 3 
per cent of radio-identified SMGs are broad-line QSOs (Chapman 
et al. 2005), suggesting that this source is not representative of the 
SMC population as a whole. SMM J04135-I-10277 is therefore not 
included in our sample of CO-detected SMGs. 

The final SMC sample consists of 12 CO detections and 6 
non-detections, their observational properties are listed in Table 2. 
The CO line profiles of the 12 detected sources are shown in Pig.Q 
In the following sections we shall use the CO observations to inves¬ 
tigate the physical properties of the SMG population, see Tabled 
In the case where a SMG is known to be gravitationally-lensed, we 
have corrected the luminosities, gas masses, and linear sizes, etc.. 
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Table 2. SMGs observed in CO to date, including the submm-bright ERO J16450+4626 which although not selected in the submm is thought to be similar to the SMG population. 
The top section of the table lists the five new CO detections presented here and the three previously published detections from our survey {Neri et al. 2003); the 6 non-detected 
SMGs are listed in the middle section, while the sources from the literature are given in the bottom section of the table. 


Source name 

Transition 

CO position 

Optical/NIR position 

Radio position 

2^spec 

•2c O 


^co 

Ref. 



<^J2000 <512000 

Q:J2000 5j2000 

Q:j2000 <5j2000 








(h ms) (° ' ") 

(h ms) (“ ' ") 

(h ms) (“ ' ") 



(kms“^) 

(Jy kms“^) 



Detections 


SMM 302396-0134“’^ 

(2-1) 

02 

39 

56.59 

-01 

34 

26.6 

02 39 56.51 

-01 34 25.66 

02 

39 56.30 

-01 

34 

30.8 

1.062 

±0.001'' 

1.062 ± 

0.002 

780 + 60 

3.4+ 0.3 

[1] 

SMM J13120+4242 

(4-3) 

13 

12 

01.20 

+42 

42 

08.8 

13 12 01.23 

+42 42 08.2 

13 

12 01.17 

+42 

42 

08.1 

3.405 

±0.001“^ 

3.408 ± 

0.002 

530 + 50 

1.7 + 0.3 

[1] 

SMMJ16366+4105 

(3-2) 

16 

36 

58.23 

+41 

05 

23.7 

16 36 58.21 

+41 05 23.9 

16 

36 58.19 

+41 

05 

23.8 

2.454 

±0.001“^ 

2.450 + 

0.002 

870 + 80 

1.8+ 0.3 

[1] 

SMMJ16371+4053 

(3-2) 

16 

37 

06.50 

+40 

53 

13.8 

16 37 06.54 

+40 53 14.1 

16 

37 06.51 

+40 

53 

13.8 

2.374 

±0.001'^ 

2.380 + 

0.004 

830 + 130 

1.0+ 0.2 

[1] 

SMM J22174+0015 

(3-2) 

22 

17 

35.20 

+00 

15 

37.6 

22 17 34.95 

+00 15 33.2 

22 

17 35.15 

+00 

15 

37.2 

3.098 

±0.002“^ 

3.099 + 

0.004 

780 + 100 

0.8+ 0.2 

[1] 

SMM 304431+0210“’'' 

(3-2) 

04 

43 

07.25 

+02 

10 

23.3 

04 43 07.25 

+02 10 24.4 

04 

43 07.25 

+02 

10 

24.4 

2.5092 

± 0.0008' 

2.5094 + 

0.0002 

350 + 60 

1.4+ 0.2 

[2] 


(7-6) 



















< 0.8 

[2] 

SMM 309431+4700“ 

(4-3) 

09 

43 

03.74 

+47 

00 

15.3 

09 43 03.70-'’ 

+47 00 15.1-'^ 

09 43 03.7 

+47 

00 

15.1 

3.349 

±0.001“^ 

3.3460 + 

0.0001 

420 + 50 

1.1 + 0.1 

[2] 


(9-8) 



















< 1.0 

[2] 

SMM 316368+4057 

(3-2) 

16 

36 

50.43 

+40 

57 

34.7 

16 36 50.40 

+40 57 34.2 

16 

36 50.43 

+40 

57 

34.5 

2.380 

±0.002“^ 

2.3853 + 

0.0014 

840+ 110 

2.3+ 0.2 

[2] 


(7-6) 

16 

36 

50.41 

+40 

57 

34.3 










2.383 + 

0.002 


1.1+ 0.2 

[2] 

Non-detections 

SMM 310523+5722 

(3-2) 







10 52 30.75 

+57 22 09.4 

10 

52 30.73 

+57 

22 

09.5 

2.611 

±0.001“^ 

2.5901 - 

2.6119 


< 0.6® 

[1] 

SMM 310524+5724 

(3-2) 







10 52 38.34 

+57 24 36.0 

10 

52 38.30 

+57 

24 

35.8 

3.036 

±0.001“^ 

3.0593 - 

3.0318 


< 1.0® 

[1] 

SMM 312360+6210 

(3-2) 







12 36 00.12 

+62 10 47.7 

12 

36 00.20 

+62 

10 

47.0 

1.994 

±0.001“^ 

1.9865 - 

2.0015 


< 1.3® 

[1] 

SMM 313123+4239 

(3-2) 







13 12 32.31 

+42 39 49.5 

13 

12 32.30 

+42 

39 

50.0 

2.320 

±0.001“^ 

2.3115 - 

2.3300 


< 1.4® 

[1] 

SMM 316363+4055 

(3-2) 







16 36 31.47 

+40 55 46.6 

16 

36 31.47 

+40 

55 

46.9 

2.283 

±0.001“^ 

2.2676 - 

2.2919 


< 0.6® 

[1] 

SMM 316363+4056 

(2-1) 







16 36 39.16 

+40 56 35.9 

16 

36 39.01 

+40 

56 

35.9 

1.495 

±0.001“^ 

1.4802 - 

1.5001 


< 1.4® 

[1] 

Literature sources 

SMM J02399-0136“'*’ 

(3-2) 

02 

39 

51.89 

-01 

35 

58.9 

02 39 51.88 

-01 35 58.0 

02 

39 51.88 

-01 

35 

58.0 

2.803 

± 0.003“^ 

2.808 + 

0.002 

710 + 80 

3.0+ 0.4 

[3] 


(3-2) 

02 

39 

51.87 

-01 

35 

58.8 

02 39 51.85 

-01 35 58.2 








2.8076 + 

0.0002 

> 1100 

3.1+ 0.4 

[4] 

SMMJ14011+0252“-*’ 

(3-2) 

14 

01 

04.92 

+02 

52 

25.6 

14 01 04.95 

+02 52 24.0 

14 

01 04.96 

+02 

52 

23.5 

2.562 

±0.002“^ 

2.5653 + 

0.0003 

200 + 40 

2.4+ 0.3 

[5] 


(3-2) 

14 

01 

04.93 

+02 

52 

24.1 










2.5652 + 

0.0001 

190 + 11 

2.8+ 0.3 

[6] 


(7-6) 

14 

01 

04.92 

+02 

52 

23.8 










2.5651 + 

0.0002 

170 + 30 

3.2+ 0.5 

[6] 

SMMJ16359+6612“'*’ 

(3-2) 

16 

35 

54.15 

+66 12 24 

16 35 54.19 

+66 12 24.9 






2.5165 

± 0.0015' 

2.5168 + 

0.0003 

500 + 100 

3.5+ 0.1 

[7] 


(3-2) 

16 

35 

54.10 

+66 12 23.8 










2.5174 + 

0.0002 

500 + 100 

2.50 + 0.12 

[8] 

ERO J16450+4626*’ 

(2-1) 

16 

45 

02.26 

+46 26 26.5 

16 45 02.36 

+46 26 25.5 






1.443 

±0.001' 

1.439 + 

0.001 

400 + 20 

1.40 + 0.10 

[8] 


(5-4) 
















1.440 + 

0.001 

380 + 20 

1.35 + 0.10 

[9] 


(1-0) 



















0.6+ 0.1 

[10] 


“ The source is lensed, see Table 3 for magnifications /xl- 
^ The CO emission is extended. 

^ The redshift is derived from the OIIA3727 line {Soucail et al. 1999). 

^ The redshift is derived from Lya (Ivison et al. 2000; Ledlow et al. 2002; Small et al. 2003; Chapman et al. 2005). 

® The redshift is derived from Ho (Dey et al. 1999; Prayer et al. 2003; Kneib et al. 2004; Swinbank et al. 2004) . 

^ Position of H7, see Neri et al. (2003) for details. 

^ The 3 <t upper limits have been calculated over an assumed velocity range of 500kms~^. 

fl] This work; [21 Neri et al. (2003); [3] Prayer et al. (1998); [4] Genzel et al. (2003); [5] Prayer et al. (1999); [6l Downes & Solomon (2003); [7l Sheth et al. (2004); [8l Kneib et al. 
(2004); [9] Andreani et al. (2000); [10] Greve et al. (2003). 
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Figure 3. The histograms show the distribution of velocity differences de¬ 
rived from the redshifts from CO line observations and restframe-UV spec¬ 
troscopy. We show the distribution for the eight SMGs for which spec¬ 
troscopic redshifts have been derived from the Lya line (grey shaded his¬ 
togram) and for the 15 QSOs in Table l4l (two of the QSOs have velocity- 
offsets ^ —lOOOkms”^). The hashed histogram compares the redshifts 
between the CO and Hck emission line for the 4 SMGs with redshifts de¬ 
rived from Ha (3) or [Oil] A3727 A (1). The bin size is 150 km s“^ which 
also reflects the typical error in the velocity offsets. 


using the gravitational lensing magnification factors listed in Table 

m 

4.2 Comparison of optical and CO redshifts 

In general, CO lines provide excellent systemic redshifts for galax¬ 
ies, by tracing the extended molecular gas distribution rather than 
the ionised gas produced by shocked outflows or accretion onto 
an AGN. The latter is typically traced by high-ionisation, e.g. 
broad ClvA1549 lines, which have been shown to be systemati¬ 
cally blue-shifted with respect to the systemic redshift in quasars 
(e.g. Richards et al. 2002), and by the Lya line, which in the ma¬ 
jority of high-redshift CO detections is blue-shifted with respect to 
the CO redshift due to outflows and dust obscuration (e.g. Mainline 
et al. 2004). The purpose of this section is to determine if similar 
systematic offsets between the CO and Lya redshifts exist within 
the SMG population. We have measurements of the Lya redshifts 
for eight of the CO-detected SMGs, see Table 2. 

A large fraction of CO line profiles shown in Fig.Qare double 
peaked or asymmetric, making it difficult to accurately determine 
the CO redshift by fitting Gaussian profiles to the spectra. Instead, 
the CO redshift of a source was determined by computing its flux- 
weighted redshift, i.e. Zco = X/ S ^(^)' error on the 

redshift is given by A^co = X/ ~ ^co)^/ X] The CO 

redshifts computed in this way are listed in Table 2. 

In Fig. 13 we have plotted the distribution of velocity offsets 
corresponding to the differences between the CO and Lya redshifts 
for the SMG sample. Negative velocity-offsets correspond to blue- 
shifted UV emission with respect to the systemic CO redshift. Also 
shown are the velocity offsets for the four SMGs with Ha redshifts. 

The velocity offsets between the Lya and CO redshifts of 


SMGs show a tendency for blue-shifted offsets (five and two SMGs 
have blue- and red-shifted offsets, respectively, with one source 
consistent with no velocity offset within the errors, ~ 150 km s“^), 
although this needs to be confirmed once a larger sample of SMGs 
is studied. The 15 QSOs detected in CO to date show a similar 
tentative excess of negative CO-Lya velocity offsets (Fig.|3, al¬ 
though in this case the the shift is due to strong absorption of the 
line - most likely from gas on small scales near the nucleus. In 
the case of SMGs, the offsets are likely to come from absorption 
from a wind on larger kpc-scales - where we only see the line scat¬ 
tered from the far side of the wind/shell due to the absorption of the 
blueshifted component. 

It is clear from Fig. |3 that in the majority of the SMGs, sig¬ 
nificant velocity offsets exist between the Lya and CO emission. 
Such large offsets illustrate how easy it would be to miss a CO line 
based on blind tuning to the restframe-UV redshift, and highlight 
the need for future submm/mm correlators with large instantaneous 
bandwidths. 

Realising the potential danger of not detecting some of the 
SMGs because of large offsets between the CO and restframe-UV 
redshifts, we have undertaken a programme to target SMGs us¬ 
ing the near-infrared spectrographs NIRSPEC on Keck-II, OHS 
on Subaru and ISAAC on the VLT, to allow us to measure red¬ 
shifts from restframe optical emission lines, which should pro¬ 
vide more reliable estimates of the systemic redshifts (e.g. Simp¬ 
son et al. 2004; Swinbank et al. 2004). The usefulness of this ap¬ 
proach is demonstrated by the hashed histogram in Fig. |3 which 
shows that for the four SMGs with redshifts derived from Ha (or 
[Oil] A3727 A in the case of SMM J02396-0134 - Soucail et al. 
1999), the velocity offsets amount to no more than ±75kms“^ . 


4.3 CO Luminosities and gas masses 

The CO line luminosities, L'co, of the individual SMGs were 
derived from their velocity-integrated line flux densities follow¬ 
ing Solomon et al. (1997), and corrected for gravitational lens¬ 
ing if necessary (see Table |3. For the eight sources detected 
in CO(3—2) we find a median luminosity of (Lco( 3 - 2 )) ~ 

(3.8 ± 2.3) X Kkms“^ pc^. Two sources were detected in 
CO(4-3) (SMM J09431-b4700 and SMM J13120-b4242) and two 
in CO(2-l) (SMMJ02399-0134 and FRO J16450-f4626). Note 
that for HR10/EROJ16450-f4626 we have used the CO(2-l) 
spectrum (Andreani et al. 2000) rather than the CO(l—0) measure¬ 
ment by Greve et al. (2003) which lacks velocity information. 

Assuming intrinsic velocity/area-averaged brightness temper¬ 
ature line ratios corresponding to an optically thick, thermalised 
gas, i.e. r 32 = Tb(3-2)/Tb{l-0) = L'(3-2)/L'(l-0) = 1, and 
similarly r 43 = 1 and r 2 i = 1, we derive a median CO(l—0) lu¬ 
minosity of = (3.8±2.0) X lO*^® Kkms”^ pc^, where 

we have included all sources, except for the non-detections. If we 
assume more realistic values of r 32 = 0.64 (Devereux et al. 1994), 
r 43 = 0.45 (Papadopoulos et al. 2000), and r 2 i = 0.9 (e.g. Braine 
& Combes 1992; Aalto et al. 1995) as derived from ISM studies in 
local starburst galaxies, we instead find a median CO(l—0) lumi¬ 
nosity of (Lco(i_o)) = (5-9 T 3-6) x 10^^° Kkms“^ pc^. 

The observed CO luminosities are converted into molecu¬ 
lar gas masses using M(H 2 ) = A'coL(,o(]^_g), where Xco = 
0.8 (Kkms“^ pc^)~^Mq is the conversion factor appropriate for 
UV-intense environments, as derived from observations of lo¬ 
cal ULIRGs (Downes & Solomon 1998). The molecular gas 
masses derived for each source are given in Table |3 for assumed 
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(2-l)/(l-0), (3-2)/(l-0) and (4-3)/(l-0) line ratios of unity. 
The distribution of SMG molecular gas masses is shown in Fig.|3 
from which we derive a median value (M(H 2 )) = (3.0 ± 1.6) x 
10^° Mq (not including the upper limits from the non-detections). 
Given the modest size of our sample we are not in a position to 
separate the observed scatter in the distribution of molecular gas 
masses into components due to varying inclination angles, differ¬ 
ences in excitation conditions, and underlying scatter in the mass 
distribution of SMGs. 

Since these molecular gas mass estimates are based on obser¬ 
vations of high-J CO transitions with J > 2 and made assum¬ 
ing a thermalised, optically thick gas, it is possible that signifi¬ 
cant amounts of cold, possibly sub-thermal, molecular gas, could 
be present, but only detectable in lower CO transitions (e.g. Pa- 
padopolous et al. 2001; Papadopoulous & Ivison 2002; Greve et 
al. 2003). Furthermore, if metal-poor gas is present, the gas mass 
could be higher. Finally, CO is primarily a tracer of material in the 
diffuse ISM and it is possible that a significant amount of dense and 
clumpy gas (<; 10"^ cm“^) is missed by our observations (e.g. Car- 
illi et al. 2004). The quoted H 2 masses in Tablel^should therefore 
be considered as lower limits to the total amount of molecular gas. 

Our estimates of the molecular gas mass include a non- 
negligible fraction of Helium which is accounted for in the adopted 
value of the conversion factor (Downes & Solomon 1998). The 
amount of neutral gas (Hi), however, is particular difficult to de¬ 
termine and significant uncertainty is associated with M(Hl)- 
estimates at high redshift (e.g. De Breuck et al. 2003a). Estimates of 
the H 2 -to-Hl mass ratio in infrared-luminous /RAS galaxies range 
from ~ 0.5 in systems with Lfir ~ Lq (Andreani et al. 

1995) to ~ 4 for Lfir ~ 10*^^ Lq systems (Mirabel & Sanders 
1989). The latter value may apply to the extremely far-infrared lu¬ 
minous SMGs. However, most of the neutral gas in local galaxies 
is on very large scales (> 10 kpc) and thus may not be relevant for 
our estimates. 

Finally, we have looked for evidence of a correlation be¬ 
tween molecular gas content and redshift. Such a trend would be 
indicative of evolution within the redshift range spanned by the 
SMG population. If we only consider the 13 sources observed in 
CO as part of our statistically-complete survey (i.e. not includ¬ 
ing SMM J02396—0134) and split the sample into low- and high- 
redshift halves at the median redshift of the sample, < 2 >= 2.4, 
we find that the CO detection rate in the high-redshift (2 > 2.4) 
subset is 5/7 (71 per cent) compared to just 2/6 (33 per cent) for 
sources at 2 < 2.4. The higher detection rate at 2 > 2.4 could 
be indicative of evolution over the redshift range 2 ~ 1-3 in the 
gas masses associated with the most luminous starbursts. Using a 
Cox/Hazard survival analysis, we find that the formal likelihood of 
a variation in detection rate with redshift in the 13 sources observed 
by our survey is P = 0.995, suggestive of a real trend. If we in¬ 
stead include all the sources listed in Table 2 the trend weakens (the 
detection rates become 4/8 and 8/10 for the low- and high-redshift 
bins, respectively). This is unsurprising since the bias against pub¬ 
lishing non-detections in the literature means that this comparison 
is less reliable than that for our statistically-complete sample. 

4.4 Line widths and dynamical masses 

A striking feature of the SMG population is their typically broad 
CO line profiles. Tabled The median observed FWHM of the sam¬ 
ple is 780 ± 320 kms”*^ and the median FWZI is 850 hms”*^. 
These values do not take into account any geometrical effects, such 
as inclination angle, which could effect the individual estimated 


line widths by more than a factor of two. For example, in the case 
of SMM J14011-1-0252 it is plausible that the relatively small line 
width is due to the fact that the angular momentum vector of the 
system is closely aligned to our line-of-sight (Tecza et al. 2004). 
Hence the observed FWHM are firm lower-limits to the true orbital 
velocities and assuming our sample has random orbital inclinations, 
then the average correction should be sec 7r/4 = 1.4 x. 

In order to use our observed CO line widths to derive the dy¬ 
namical mass of the SMGs we must know the spatial extent of the 
CO. However, so far CO emission has been reliably resolved in 
only one SMG: SMM 102399-0136 (Genzel et al. 2003). While a 
number of SMGs show tentative evidence of extended CO emis¬ 
sion (Ivison et al. 2001; Neri et al. 2003; Greve et al. 2003), their 
low signal-to-noise ratios prevent a robust determinations of the 
source size. Instead, we have assumed a conservative source di¬ 
ameter of ~ 0 . 5 ”, which corresponds to a median linear diameter 
of 3.7kpc (correcting for gravitational lensing where relevant), see 
Tabled This is a factor of two smaller than the source sizes esti¬ 
mated by Neri et al. (2003), but consistent with characteristic size 
estimates from recent high-resolution radio observations of SMGs 
with MERLIN by Chapman et al. (2004), and with results from 
high-resolution PdBI CO observations (Tacconi et al. 2005) of the 
three SMGs detected in CO by Neri et al. (2003). 

Another notable feature of our sample is the high fraction 
of SMGs which show evidence of double-peaked CO line pro¬ 
files. Four sources (SMM 102396-0134, SMM 102399-0136, 
SMM 116359-f 6612 and SMM 1163664-4105) show unam¬ 
biguous double-peaked profiles, and a further two sources 
(SMM104431-f0210 and SMM 116368-44057) are better fit by 
a double Gaussian than a single Gaussian, as measured by the 
reduced of the fit. In the case of SMM 1163684-4057 further 
evidence of multiple CO peaks comes from Ha IFU observations 
(Swinbank et al. in preparation) which show emission components 
at velocity-offets similar to those we infer from the double 
Gaussian fit. Thus, at least 4/12 (33 per cent), and possibly as 
many as 6/12 (50 per cent), of the detected sample show evidence 
of having more than one CO-emitting component in their spectra. 
Comparing the occurrence of double-peaked line profiles in 
sources as a function of redshift within our sample, we see that 
the double-peaked sources appear to lie at the lower-redshift end 
of our sample, 5/7 at 2 < 2.5, compared to 1/5 for the SMGs at 
2 > 2.5. A possible explanation for this tentative trend is discussed 
later. 

Such multi-peaked line profiles are a tell-tale sign of orbital 
motion under the influence of gravity, and can usually be attributed 
to either a disk or a merger. To determine the likely structure of 
the gas reservoir, we use our median gas mass (3 x Mq ) and 
adopt a typical diameter of ~ 4 kpc for the gas reservoir, to derive a 
mass surface density of (E ~ 2.4 x 10® Mq pc“^) for this struc¬ 
ture. This is extremely high and implies that if the gas is present 
in a disk then this will have a Toomre parameter Q < 1®, indicat¬ 
ing it will be unstable to bar formation. The disk would then col¬ 
lapse on a timescale comparable to the sound crossing time, which 
is much shorter than the expected ^ 10-100 Myr duration of the 
SMG phase, see 45.41 We conclude that the kinematics responsi¬ 
ble for the double-peaked CO lines are unlikely to arise from gas 


^ Toomre’s stability criterion says that Q = > 1 in order for a 

gaseous disk to be stable: see e.g. Binney & Tremaine (1987). For a surface 
density S ~ 2.4 X 10® Mq pc~^, Q will always be less than one for any 
realistic values of the sound speed Vs and circular frequency Q. 


© 0000 RAS, MNRAS 000, 000-000 


10 Greve et al. 


Table 3. Physical properties derived from the CO observations. 


Source 


(kpc) 

Transition 

r / b 

^CO 

(xl0l”Kkms“lpc^) 

M(H2)'>'‘= 

(Xl030 Mq ) 

AV^ 

^ ^FWHM 

(kms~^) 

m! 

dyn 

(xIO^^Mq ) 

SMM J02396-0134 

2.5 

3.3 

(2-1) 

2.1 ± 0.2 

1.7 ± 0.2 

780 ± 60 

1.2 

SMMJ02399-0136 

2.5 

3.2 

(3-2) 

4.8 ± 0.8 

3.8 ± 0.6 

1360 ± 50 

3.5 

SMM J04431-I-0210 

4.4 

1.9 

(3-2) 

1.1 ± 0.2 

0.9 ± 0.2 

350 ± 60 

0.1 

SMMJ09431-1-4700 

1.2 

6.3 

(4-3) 

2.7 ± 0.3 

2.2 ± 0.2 

420 ± 50 

0.7 

SMMJ10523+5722 

1.0 

8.1 

(3-2) 

< 2.1 

< 1.4 



SMMJ10524-I-5724 

1.0 

7.8 

(3-2) 

< 4.6 

< 3.0 



SMMJ12360-1-6210 

1.0 

8.5 

(3-2) 

< 2.9 

< 1.6 



SMMJ13120-1-4242 

1.0 

7.5 

(4-3) 

5.3 ± 0.9 

4.2 ± 0.7 

530 ± 50 

1.2 

SMMJ13123-I-4239 

1.0 

8.3 

(3-2) 

< 4.1 

< 2.4 



SMMJ14011-(-0252'f 

5.0 

1.6 

(3-2) 

1.9 ± 0.2 

1.5 ± 0.2 

190 ± 11 

0.03 

SMMJ16368-1-4057 

1.0 

8.3 

(3-2) 

7.0 ± 0.6 

5.6 ± 0.5 

840 ± 110 

3.5 

SMMJ16359-i-6612 

22 

0.4 

(3-2) 

0.4 ± 0.2 

0.3 ± 0.2 

500 ± 100 

0.06 

SMM 116363+4055 

1.0 

8.3 

(3-2) 

< 1.7 

< 0.9 



SMMJ16363+4056 

1.0 

8.5 

(2-1) 

< 4.1 

< 1.8 



SMMJ16366+4105 

1.0 

8.2 

(3-2) 

5.7 ± 1.0 

4.6 ± 0.7 

870 ± 80 

3.7 

SMMJ16371+4053 

1.0 

8.3 

(3-2) 

3.0 ± 0.9 

2.4 ± 0.7 

830 ± 130 

3.4 

EROJ16450+4626 

1.0 

8.5 

(2-1) 

3.8 ± 0.3 

3.0 ± 0.2 

400 ± 20 

0.8 

SMMJ22174+0015 

1.0 

7.8 

(3-2) 

3.8 ± 1.0 

3.0 ± 0.7 

780 ± 100 

2.8 

Median^ 


7.5 ± 3.1 


3.8 ± 2.0 

3.0 ± 1.6 

780 ± 320 

1.2 ± 1.5 


“ Assuming equal flux magnification and linear magnification. 

^ D(l^^) is the linear distance in kpc corresponding to l" at the given redshift and corrected for the lensing amplification /j,l. 
^ Derived assuming Xco=0-8(Kkms~^ pc^)~^Mq . 

The line widths are derived from a single Gaussian fit to the line profile. 

® Calculated adopting a source diameter of 0.5. 

^ The amplification factor for this source was recently revised from 2.5 to ~ 5 (see discussion in Swinbank et al. 2004). 

^ The non-detections (upper limits) are not included in the median. 



m(hO [io'° m©] 

Figure 4. The distiibutions of molecular gas masses of the SMG sample 
assuming a conversion factor of Xco=0.8(Kkms“^ pc^)“^M 0 and line 
ratios r 32 = r 43 = 1. Non-detections are included as the shaded part of 
the histogram. The median gas mass of the sample, not including the non¬ 
detections, is (M(H 2 )) = (3.0 ± 1.6) X Mq (dotted vertical line). 


distributed in a stable disk, instead a more likely scenario is that 
they reflect a merger of two gas-rich components or from a disk 
collapsing under gravitational instability. 

With almost all of the SMGs being spatially unresolved in 
CO we are unable to determine whether SMGs typically contain 
gaseous disks in their centres or if they are mergers. However, given 
their large gas mass densities, as well as the high fraction of double- 
peaked profiles we feel that SMGs are more likely to be mergers. 


Additional evidence for a merger origin comes from the detection 
of close companions in CO in some of our fields (Neri et al. 2003), 
as these these are likely to be interacting or merging given the small 
spatial and velocity differences. Thus, dynamical masses were cal¬ 
culated using the ’merger’ formula given in Neri et al. 2003 (see 
also Genzel et al. 2003), and reported in Table|3 Neri et al. (2003) 
assumed a stable disk model which gives dynamical masses lower 
by a factor of two; but since they adopted a characteristic source 
sizes twice as large as ours, the resulting dynamical mass estimates 
coincide. The median SMG dynamical mass within R ^ 2kpc is 
(Mdyn) — (1.2 db 1.5) X 10^^ Mq , where we have applied an 
average correction of 1.4x, assuming random inclinations for the 
sample overall. We estimate a median molecular gas-to-dynamical 
mass fraction in SMGs of (Mgas/Mdyn) ~ 0.3 (~ 0.5 if we as¬ 
sume a disk model). This indicates that the contribution of the gas 
reservoir to the dynamics of the central regions of typical SMGs 
is substantial, and that the gas has already concentrated into the 
central regions of the galaxy by some dissipative process. It should 
be noted, however, that significant uncertainties are associated with 
the dynamical mass estimates, especially if the system is not in dy¬ 
namical equilibrium, or if the observed line is not a good tracer of 
the dynamical state of the galaxy (also see the discussion in Neri et 
al. 2003). 


5 COMPARISON WITH OTHER POPULATIONS 

Unlike optical/near-infrared and X-ray observations, CO observa¬ 
tions are relatively unaffected by either dust extinction or AGN ac¬ 
tivity, and a comparison between the CO properties of SMGs and 
that of other galaxy populations could therefore potentially provide 
us with a clean way of comparing the bulk properties of SMGs with 
other galaxy populations at both high and low redshifts; although, 
by virtue of the submm selection technique such a comparison is of 
course biased towards dust-enshrouded sources. 
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Figure 5. The relationship of logjQ(LpQ) with redshift for the (U)LIRGs 
samples of Sanders et al. (1991), Solomon et al. (1997), and Yao et 
al. (2003), and for the sample of high-redshift (sub)mm-selected sources 
presented in this paper. The CO luminosities have been corrected for 
gravitational amplification where necessary. Also shown are the (lensing- 
conected) CO luminosities for HzRGs, QSOs and LBGs listed in Table|4] 
with the exception of Q 0957+561 and MG 0414+0534 which have un¬ 
known gravitational magnification factors. The CO luminosities for the 
(U)LlRGs are all based on the CO(l—0) line, whereas the luminosities 
of the high-redshift sources are entirely derived from CO J = J + 1 —► J 
transitions with J + 1 > 2 (Table |4}. Upper limits on the CO luminosi¬ 
ties for the 7 SMGs not detected in CO are represented by downward 
pointing arrows. The hashed area denotes the region where sources are 
precluded from detection due to an integrated CO(3 — 2) flux limit of 
Ico = 0.3Jykms“^. Of our sample, only SMMJ16359+6612 with its 
very large amplification factor (p, = 22) is detectable below this limit. 


We list in Table^a// non-SMG CO sources detected at 2 > 1. 
The bulk of these are extreme AGN such as submm-bright, high- 
redshift radio galaxies (HzRGs) or QSOs (e.g. Omont et al. 1996; 
Papadopoulos et al. 2000; Cox et al. 2002). 

Due to their similar properties, in particular their large far- 
infrared luminosities, distorted morphologies and optical/near- 
infrared colours (e.g. Ivison et al. 2002; Small et al. 2004), SMGs 
are commonly thought to be high-redshift analogues of local 
ULIRGs. To test this claim we now compare the CO properties 
of SMGs with ULIRGs. To this end we have used the samples of 
Sanders et al. (1991; SA91) and Solomon et al. (1997; S097) which 
consist of 48 and 37 local ULIRGs (and less luminous LIRGs) 
in the redshift range 2 = 0.03-0.27, respectively, and the sam¬ 
ple of 60 (U)LIRGs from Yao et al. (2003; Y03), selected from 
the SCUBA Local Universe Galaxy Survey (SLUGS, Dunne et al. 
2000 ). 


5.1 CO luminosities and gas masses 

In Fig. 13 we plot the CO luminosities of the SMG sample as a 
function of redshift along with the CO(l—0) luminosities of the 
(U)LIRG samples of SA91, S097, and Y03. From Fig. 0 it is 


clear that while the most luminous ULIRGs have CO luminosi¬ 
ties comparable to those of the faintest detected SMGs, SMGs are 
generally more CO luminous than the local (U)LIRGs. In 1|4.3l we 
found the median CO line luminosity of the SMG sample to be 
(Leo) = (3.8 ± 2.0) X 10^°Kkms“’^ pc^, which is almost a factor 
of four greater than the average CO luminosity of the ULIRGs from 
S097. Furthermore, while the CO luminosities of the (U)LIRGs 
are all based on the lowest and least excitation-biased CO(l—0) 
rotational line, all of the SMGs are detected in higher transitions. 
Hence, the CO luminosities based on the high-J lines for SMGs 
should be regarded as lower limits on the CO(l—0) luminosities, 
unless the lines are fully thermalised at high temperatures and op¬ 
tically thick. 

The (U)LIRGs lie on a well-defined locus in Fig.|3 which is 
extended to the highest redshifts and CO luminosities by the SMG 
sample. The lower boundary of this trend is determined by the sen¬ 
sitivity limits of the CO observations, while the upper bound is a 
real astrophysical limit. The steady decline in the CO luminosity of 
SMGs to ULIRGs and LIRGs seen in Fig.|3is likely to reflect the 
evolution in the molecular gas-content of the most luminous star- 
burst galaxies as a function of redshift, despite the expected increas¬ 
ing mass of a typical galaxy with cosmic time. A possible caveat 
is the assumption of a constant conversion factor in (U)LIRGs and 
SMGs. In (U)LIRGs the Xco conversion factor is about 4.5 times 
lower than in normal spiral galaxies, and it is possible that for the 
yet more luminous SMGs it is even lower. 

Also shown in Fig. |3 are the CO luminosities of all HzRGs 
and QSOs detected to date (Table|3, as well as the only optically- 
selected Lyman break galaxy (LBG), MS 1512—cB58 (z = 2.73), 
detected in CO (Baker et al. 2004). Their large apparent scatter 
Is due to many of these objects being gravitationally-lensed, in¬ 
creasing the effective sensitivity limits of their observations beyond 
the capabilities of current instruments in blank fields. SMGs are 
seen to have CO luminosities comparable to the most luminous 
HzRGs (e.g. Papadopoulos et al. 2000; De Breuck et al. 2003a, 
2003b) and QSOs (e.g. Carilli et al. 2002a; Walter et al. 2003). In 
contrast, the CO luminosity of the typical-luminosity LBG cB58, 
4.2 X 10®Kkm s“^ pc^, is nearly two orders of magnitude (^ 90x) 
less than the median value for SMGs, despite its similar redshift. 
With only one LBG detected in CO so far it is not possible to 
make any meaningful conclusions about the molecular gas content 
of this population, except that the low CO luminosity of cB58 com¬ 
pared to that of SMGs is consistent with the faintness of LBGs at 
submm wavelengths and their relatively low (compared to SMGs) 
star-formation rates as inferred from optical spectroscopy (Blain et 
al. 1999; Adelberger & Steidel 2000). 


5.2 CO line widths and dynamical masses 

In Il4.4l we found that ^ 33 per cent of the CO-detected SMGs have 
double-peaked CO profiles. This is comparable to the fraction of 
double-peaked CO spectra in local ULIRGs (S097). The median 
FWHM CO linewidth for the SMG sample is 780 ± 320kms“^ , 
about 3 and 4 times larger than the averages for the SA91 and S097 
ULIRG samples, respectively. Since the dynamical mass depends 
on the square of the velocity dispersion (Mdyn ^ Ro^), this would 
naively suggest that the dynamical masses of the SMGs are ~ 9- 
16 times larger than that of ULIRGs (if the gas kinematics sample 
similar radii in the galaxies). From detailed interferometric stud¬ 
ies of the circumnuclear gas in local ULIRGs Downes & Solomon 
(1998) find a median dynamical mass of ~ 6 x 10® Mq within 
R ^ 0.6 kpc - about 20x smaller than the mass enclosed within 


© 0000 RAS, MNRAS 000, 000-000 



12 Greve et al. 


Table 4. List of published high-redshift CO detections of LBGs, QSOs and HzRGs. The sources are sorted according to their redshifts. This table, together 
with Table 2summarises the complete list of CO detections of sources at 2 : > 1 to date. The velocity-integrated line fluxes, Jco, have not been corrected for 
gravitational amplification. 


Source 

Type 

Transition 

•2^spec 

2CO 

AVpWHM 

(kms~^) 

^co 

(Jykms“^) 

Ref. 

MS 1512-cB58“ 

LBG 

(3-2) 

2.727 

2.7265 ± 0.0004 

174 ± 43 

0.37 ± 0.08 

[1] 

Q0957-|-561“ 

QSO 

(2-1) 

1.413 

1.414 

440 

1.2 

[2],[3] 

IRASF102I4+4724“ 

QSO 

(3-2) 

2.286 

2.2867 ± 0.0003 

250 

21 

[4],[5] 



(3-2) 


2.2855 ± 0.0003 

230 ± 30 

4.1 ± 0.9 

[6] 



(3-2) 


2.2854 ± 0.0001 

220 ± 30 

4.2 ± 0.8 

[7] 



(6-5) 


2.2857 ± 0.0003 

240 ± 30 

9.4 ± 2.0 

[6] 

H 1413+117“ 

QSO 

(3-2) 

2.5582 ± 0.0003 

2.558 

326 

8.1 

[8] 



(3-2) 



352 ± 81 

14.4 ± 4.4 

[9] 



(3-2) 


2.5579 

362 ± 23 

9.9 ± 0.6 

[10] 



(3-2) 


2.55784 ± 0.00003 

416 ± 6 

13.2 ± 0.2 

[11] 



(4-3) 


2.5579 

375 ± 16 

21.1 ± 0.8 

[10] 



(5-4) 


2.5579 

398 ± 25 

24.0 ± 1.4 

[10] 



(7-6) 


2.5579 

376 

47.3 ± 2.2 

[10] 

VCVJ1409+5628 

QSO 

(3-2) 

2.562 

2.585 ± 0.001 

370 ± 60 

2.4 ± 0.7 

[12] 



(3-2) 


2.5832 ± 0.0001 

311 ± 28 

2.3 ± 0.2 

[13] 



(7-6) 




4.1 ± 1.0 

[13] 

MG0414+0534“ 

QSO 

(3-2) 

2.639 ± 0.002 

2.639 

580 

2.6 

[14] 

LBQS 1230+1627B 

QSO 

(3-2) 

2.735 ± 0.005 

2.741 ± 0.002 


0.80 ± 0.26 

[15] 

RX10911+0551“ 

QSO 

(3-2) 

2.800 

2.796 ± 0.001 

350 ± 60 

2.9 ± 1.1 

[12] 

SMM104135+1027“ 

QSO 

(3-2) 

2.837 ± 0.003 

2.846 ± 0.002 

340 ± 120 

5.4 ± 1.3 

[12] 

MG0751+2716“ 

QSO 

(4-3) 

3.200 ± 0.001 

3.200 

390 ± 38 

5.96 ± 0.45 

[16] 

APM 08279-1-5255“’^ 

QSO 

(1-0) 

3.87 

3.9 


0.150 ± 0.045 

[17] 



(1-0) 




0.22 ± 0.05 

[18] 



(2-1) 




1.15 ± 0.54 

[17] 



(4-3) 


3.9114 ± 0.0003 

480 ± 35 

3.7 ± 0.5 

[19] 



(9-8) 


3.9109 ± 0.0002 


9.1 ± 0.8 

[19] 

PSS12322+1944“ 

QSO 

(1-0) 

4.1108 ± 0.0005 

4.1192 ± 0.0004 

200 ± 70 

0.19 ± 0.08 

[20] 



(2-1) 




0.92 ± 0.03 

[20] 



(4-3) 


4.1199 ± 0.0008 

375 ± 41 

4.21 ± 0.40 

[21] 



(5-4) 


4.1199 ± 0.0008 

273 ± 50 

3.74 ± 0.56 

[21] 

BRI1335-0417^ 

QSO 

(2-1) 

4.398 ± 0.028 

4.4074 ± 0.0015 

420 ± 60 

0.44 ± 0.08 

[22] 



(5-4) 




2.8 ± 0.3 

[23] 

BRI0952-0115“ 

QSO 

(5-4) 

4.426 ± 0.020 

4.4337 ± 0.0006 

230 ± 30 

0.91 ± 0.11 

[15] 

BRI1202-0725^ 

QSO 

(2-1) 

4.69 



0.49 ± 0.09 

[22] 



(4-3) 



280 ± 30 

1.5 ± 0.3 

[24] 



(5-4) 


4.6932 ± 0.002 

320 ± 35 

2.40 ± 0.30 

[25] 



(5-4) 



220 ± 74 

2.7 ± 0.41 

[25] 



(7-6) 


4.6915 ± 0.001 

250 - 300 

3.1 ± 0.86 

[24] 

SDSS11148+5251 

QSO 

(3-2) 

6.43 ± 0.05 

6.418 ± 0.004 

320 

0.18 ± 0.04 

[26] 



(6-5) 


6.4187 ± 0.0006 

279 

0.73 ± 0.076 

[27] 



(7-6) 


6.4192 ± 0.0009 

279 

0.64 ± 0.088 

[27] 

53W002 

HzRG 

(3-2) 

2.390 ± 0.004 

2.394 ± 0.001 

540 ± 100 

1.5 ± 0.2 

[28] 



(3-2) 


2.3927 ± 0.0003 

420 ± 40 

1.20 ± 0.15 

[29] 

B3J2330-I-3927'’ 

HzRG 

(4-3) 

3.087 ± 0.004 

3.094 

500 

1.3 ± 0.3 

[30] 

TN10121 + 1320 

HzRG 

(4-3) 

3.516 

3.520 

700 

1.2 ± 0.4 

[31] 

6C 1909+722 

HzRG 

(4-3) 

3.5356 

3.532 

530 ± 70 

1.62 ± 0.30 

[32] 

4C41.17'’ 

HzRG 

(4-3) 

3.79786 ± 0.0024 

3.7958 + 0.0008 

1000 ± 150 

1.8 ± 0.2 

[33] 

4C 60.07’’ 

HzRG 

(1-0) 

3.788 

3.791 

> 1000 

0.24 ± 0.03 

[34] 



(4-3) 


3.791 


2.50 ± 0.43 

[32] 

TN 10924-2201 

HzRG 

(1-0) 

5.1989 ± 0.0006 

5.202 ± 0.001 

250 - 400 

0.087 ± 0.017 

[35] 



(5-4) 



200 - 300 

1.19 ± 0.27 

[35] 

“ The source is gravitationally lensed. 








^ Tlie CO emission is resolved. 








[I] Baker et al. (2004); [2] Planesas et al. (1999); [3] Krips et al. (2003); [4] Brown & Vanden Bout (1991); [5] Solomon et al. (1992a); [6] Solomon et al. (1992b); 


[7] Downes et al. (1995); [8] Barvainis et al. (1994); [9] Wilner et al. (1995); [10] Barvainis et al. (1997); [11] Weiss et al. (2003); [12] Hainline et al. (2004); 



[13] Beelen et al. (2004); [14] Bai'vainis et al. (1998); [15] Guilloteau et al. (1999); [16] Barvainis et al. (2002); [17] Papadopoulos et al. (2001); [18] Lewis et al. (2002); 
[19] Downes et al. (1999); [20] Carilli et al. (2002b); [21] Cox et al. (2002); [22] Carilli et al. (2002a); [23] Guilloteau et al. (1997); [24] Omont et al. (1996); 

[25] Ohta et al. (1996); [26] Walter et al. (2003); [27] Bertoldi et al. (2003); [28] Scoville et al. (1997); [29] Alloin et al. (2000); [30] De Breuck et al. (2003a); 

[31] De Breuck et al. (2003b); [32] Papadopoulos et al. (2000); [33] De Breuck et al. (2005); [34] Greve et al. (2004a); [35] Klamer et al. 2005. 


R < 2kpc in SMGs (see Table|3. Assuming an isothermal mass 
distribution, the SMGs will have dynamical masses within 0.6 kpc 
roughly 6x that of local ULIRGs. 

We now turn to the comparison between SMGs and other 
high-redshift galaxy populations. Fig.|^ compares the distribution 
of velocity dispersions of our SMG sample with that of all high- 
redshift QSOs and HzRGs detected in CO to date. While there is 


a substantial overlap between the SMG and QSO distributions, the 
latter is seen to peak at somewhat lower velocity dispersions than 
SMGs and we note that all QSOs have log cr < 2.4 (correspond¬ 
ing to FWHM < 590 kms“^ ). Also, no QSO detected in CO show 
evidence of a double peaked line profile (although a few are spa¬ 
tially resolved - Papadopoulos et al. 2001; Carilli et al. 2002a). A 
Kolmogorov-Smirnov test suggests that the two distributions are 
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Figure 6. a) Distiibution of CO velocity dispersions for high-redshift QSOs 
and HzRGs from Table 2 and our sample of CO-detected SMGs. In the 
cases of the two HzRGs 4C 60.07 and TN J0924—2201 we have adopted 
FWHM linewidths of 1000 and 325 km s“^ , respectively. In b), we compare 
with the measurements of the line widths of the Ha line from near-infrared 
spectroscopy of a large sample of SMGs (Swinbank et al. 2004). 


formally different at the 95 per cent confidence level. The appar¬ 
ently lower CO linewidths of QSOs is not clear, although it is pos¬ 
sible that they are either intrinsically lower-mass systems or have 
smaller gas disks. In contrast, there is no discemable difference be¬ 
tween the CO line widths of the SMGs and those HzRGs detected 
in CO. This suggest that the brightest SMGs have similar dynami¬ 
cal masses, and therefore possibly similar dark-matter halo masses, 
to HzRGs which are believed to be amongst the most massive ob¬ 
jects in the high-redshift Universe (inferred from both the presence 
of a supermassive black hole and the strong clustering of HzRGs - 
e.g. Ford et al. 1994; Kooiman et al. 1995). However, such a con¬ 
clusion is uncertain due to the small number of objects involved and 
will have to wait until a larger, more uniform sample of HzRGs has 
been observed in CO. 



lOQio {LpiR [L0]) 

Figure 7. A comparison of the CO and far-infrared luminosities 
Tfir) for SMGs and local ULIRGs. As in Fig. the an'ows denote the 
6 non-detections. The solid and dashed lines represent fits of the form 
log Lpg = o log Lfir -f /3 to all three (U)LIRGs samples and to the com¬ 
bined (U)LIRGs and SMGs samples, respectively. The dotted line is the 
best fit to all four samples but where the slope has been fixed to unity. 


In Fig.|^ we compare the CO line widths from our SMG sam¬ 
ple with measurements of the width of the Ha line from a recent 
near-infrared spectroscopic survey of 28 SMGs and high redshift, 
optically-faint, radio-selected starburst galaxies (OFRGs — see e.g. 
Chapman et al. 2005) by Swinbank et al. (2004). The two distribu¬ 
tions span a similar range of velocity dispersions, but the CO sam¬ 
ple is skewed towards higher values, resulting in a higher median 
value. We find that the average ratio between the CO and Ha line 
widths for SMGs is (fwhMco) / (fwhMhq) ~ 2. The main reason 
for this is likely to be due to the CO tracing multiple components 
in the SMG whereas the long-slit Ha spectroscopy is identifying 
only a single component due to the slit orientation or due to ex¬ 
treme obscuration in the second component. Part of the discrep¬ 
ancy may also reflect dust-obscuration within the galaxy providing 
only a partial view of the kinematics of the system, in contrast to 
the obscuration-independent measurement provided by CO. If this 
is the case, it could mean that the dynamical masses of SMGs are 
larger than suggested from their rest-frame optical spectra. 

5.3 Star-formation efficiency and the L'co-Lum relation 

While the luminosity of the CO lines is an indicator of how much 
molecular gas a galaxy contains, it does not necessarily mea¬ 
sure how efficiently this gas is being turned into stars. The star- 
formation efficiency (SFE) can be estimated from the ratio of the 
far-infrared luminosity of the system to the amount of molecular 
gas available to form stars, i.e. SFE= Lfir/M(H 2 ). 

The far-infrared luminosities of the SMGs were estimated fol¬ 
lowing Neri et al. (2003), Lfir = 1.9 x lO^^S'sso, where Lfir 
is in Lq and Ssso in mjy. This assumes a modified grey-body 
model for the far-infrared emission with a dust temperature of 
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Td = 40 K and emissivity oc The median far-infrared lu¬ 

minosity of the SMG sample is (Lfir) = (1-5 ± 0.7) x lO*^® Lq 
which is nearly an order of magnitude larger than the most lumi¬ 
nous local ULIRGs (Solomon et al. 1997) . Assuming that the hulk 
of the far-infrared luminosity is powered by star-formation (Prayer 
et al. 1998; Alexander et al. 2005), we find a median star-formation 
efficiency of (SFE) = 450 ± 170 Lq Mq for the SMG sam¬ 
ple, in agreement with the initial findings of Neri et al. (2003). This 
is somewhat higher than the typical starformation efficiency of the 
local ULIRGs studied by S097 ( {SFE) = 180 ± 160 Lq Mq-^ 
- using a conversion factor of Xco = 0.8 (Kkm s“^ pc^)“^M 0 ). 

Perhaps a more straightforward measure of the star-formation 
efficiency is the continuum-to-line luminosity ratio, Lfir/Lco> 
since it does not depend on Xco- Locally, (U)LIRGs are observed 
to follow a scaling relation between L'co and Lfir with the more far- 
infrared luminous galaxies having proportionally higher CO lumi¬ 
nosities (Rickard & Harvey 1984; Young et al. 1984, 1986; Sanders 
& Mirabel 1985). In Fig.^we have plotted the SMGs onto the Lco- 
Lfir diagram along with the (U)LIRGs from our three low-redshift 
comparison samples. The SMGs extend the general trend seen for 
the local (U)LIRGs out to far-infrared luminosities ^ 10^® Lq . 
A power-law fit to all three (U)LIRG samples yields log L'co = 
(0.62 ± 0.09) log Lfir -h (2.41 ± 1.06). Including the SMGs in 
the fit yields log Leo = (0.62 ± 0.08) log Lfir + (2.33 ± 0.93), 
i.e. virtually no change in the fit at all. The combined SMG and 
local ULIRG samples have a statistically-significant correlation, as 
shown by a Spearman’s rank-order correlation test, which yields a 
probability of P < 0.0001 that a random (uncorrelated) data set 
could result in the observed correlation coefficient (rs = 0.93). 
From Fig. 0 there is even some evidence of a Lcq-Lfir correla¬ 
tion within the SMG population itself (the probability of obtaining 
the observed correlation, rs = 0.70, by chance is P < 0.0142). 
Moreover, the slope of the L'co — Lfir correlation inferred from the 
fit is significantly less than unity. This is clearly illustrated in Fig. 
Qwhere the locus defined by the data points is seen to have a shal¬ 
lower slope than the line of equality. 

The observed slope of the correlation in Fig. Q implies that 
the Lfir/Lco ratio increases with Lfir. Several studies have already 
shown that ULIRGs have higher Lfir/Lco ratios than the less lu¬ 
minous LIRGs, which in turn have higher ratios than spiral galax¬ 
ies and GMCs (e.g. Sanders et al. 1986; Solomon & Sage 1988). 
This trend is shown in Fig.|^where Lfir/L/o has been plotted as 
a function of redshift for the LIRG, ULIRG and SMG samples. 
The median Lf,r/L/o ratios for LIRGs (Y03) and ULIRGs (S097) 
are 50 ± 30 and 160 ± 130 Lq (Kkms“^pc^)“^, respectively. 
The median far-infrared to CO luminosity ratio for the SMGs is 
360 ± 140 Lq (Kkms“^pc^)“^, which is about a factor of two 
higher than the median value for ULIRGs, and suggests that al¬ 
though some overlap exists between the two populations, in gen¬ 
eral SMGs have higher star-formation efficiencies, or a larger con¬ 
tribution to their far-infrared luminosities from a source other than 
star formation, e.g. an AGN. Although AGN are frequently found 
in SMGs using extremely deep X-ray observations, these observa¬ 
tions also suggest that the AGN are bolometrically insignificant in 
the vast majority of SMGs (Alexander et al. 2003, 2005). 

We underline that the above findings are strongly biased to¬ 
wards luminous SMGs, as they are based on a sample of pre¬ 
dominantly bright SMGs, This selection effect in conjunction with 
the CO detection limit of our survey prevent us from probing the 
lower right portion of Fig.[^ where objects with low Lfir/L/q ra¬ 
tios reside. It is possible that such objects, which are found lo¬ 
cally (see Fig. [^, will also be uncovered at high redshifts by fu- 



Figure 8. The star-formation efficiency as indicated by the observed ratio 
Lfir/L/q as a function of redshifts. Upward pointing arrows represent the 
6 SMGs not detected in CO. The hashed region represents the no-detection 
zone for sources with integrated CO(3 — 2) fluxes of Ico < 0.3 Jykms“^ 
and submm fluxes of Sssoiim < 4 mjy. 


ture deep submm and CO observations of a complete sample of 
high -2 starburst galaxies. Intriguingly, however, the Lfir/L/q ra¬ 
tio of SMM1163594-6612, which has a far-infrared and CO lu¬ 
minosity comparable to the least luminous ULIRGs in the S097 
sample, is higher than almost all the sources in that sample 
(360 Lq (Kkms“^pc^)“^). 

As a final check we repeated the above analysis using SMG 
far-infrared luminosities derived from their radio fluxes and using 
the radio-FIR correlation (Condon 1992; Yun, Reddy & Condon 
2001). There is tentative evidence that this correlation applies out 
to high redshifts (Garrett 2002; Appleton et al. 2004), and arguably 
provides a more reliable estimate of the far-infrared luminosity than 
that based on the submm flux which has significant uncertainty as¬ 
sociated with it due to the unknown dust properties of SMGs. The 
far-infrared luminosities derived using the radio information are 
generally larger, and the L/q — Lfir correlation therefore shallower 
(a — 0.59 ± 0.08). In this case, the SMG sample also exhibits a 
significant correlation, although this can partly be explained by the 
luminosity distance-squared stretching of the data points. 

5.4 Discussion 

In the previous three sub-sections it was found that SMGs have on 
average four times larger CO luminosities, 3-4 times larger line 
widths and a median Lfir /Uo ratio twice that of local ULIRGs. 
As argued in 44.41 the large line widths and molecular gas masses 
seen in SMGs are difficult to reconcile with a scenario in which the 
gas resides in a circumnuclear disk. This is in contrast to what is 
observed in local ULIRGs where the bulk of the gas is found in a 
R ~ 0.5 kpc disk or ring. In SMGs the gas is more likely to be 
distributed on scales of 2-3 kpc, as suggested by their radio mor¬ 
phologies (Chapman et al. 2004) as well as by high-resolution CO 
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observations of SMGs (Genzel et al. 2003; Tacconi et al. 2005). The 
typical gas-to-dynamical mass fraction in SMGs was estimated in 
1 14.41 to be ~ 0.3 - almost two times higher than the median gas 
mass fraction in ULIRGs (~ 0.16, Downes & Solomon 1998), as¬ 
suming that the CO-H 2 conversion factor is the same for ULIRGs 
and SMGs. The local gas mass fraction in ULIRGs is measured 
within the molecular disk {R ^ 0.5 kpc) - but when averaged over 
scales of 2-3 kpc, which are the typical scales probed by our CO 
observations of SMGs, the gas fraction is likely to be significantly 
smaller. Although, we caution that some of the above arguments 
are based on a relatively small number of sources detected at low 
signal-to-noise or have non-negligible uncertainties associated with 
them (such as the CO conversion factor), they do seem to sug¬ 
gest, together with the generally higher far-infrared luminosities 
of SMGs, that the latter are neither high-redshift replicas of local 
ULIRGs, nor simply scaled-up versions. Rather SMGs appear to 
more gas-rich and more efficiently star-forming than local ULIRGs. 

In this respect it is important to draw attention to observa¬ 
tional studies of the depletion of molecular gas with starbursts in 
ULIRGs. In a sample of more than 50 (U)LIRGs, Gao et al. (1999) 
found a correlation between the CO(l—0) luminosity and the pro¬ 
jected separation of the merging nuclei, which they took as evi¬ 
dence for molecular gas being rapidly depleted due to intense star 
formation as the merger progresses. If a similar picture applies at 
high redshift, the large CO luminosities and gas masses we find for 
SMGs would imply that these systems are extended, and caught 
in the early stages of merging. However, this might be a premature 
conclusion since Rigopoulou et al. (1999) did not find any evidence 
for a correlation between the gas mass and merger phase for local 
ULIRGs. 

If the above trends are verified by future observations, it is in¬ 
teresting to speculate what physical mechanisms could be respon¬ 
sible for the differences between local ULIRGs and SMGs? Both 
populations appear to result from mergers and strong interactions 
(e.g. Sanders et al. 1991; Chapman et al. 2003c), so what makes the 
starbursts so much more efficient in SMGs? Numerical simulations 
have compared the gas flow in major mergers between two galax¬ 
ies with strong bulges, and between two, gas-rich but bulge-less 
galaxies (Barnes & Hernquist 1996; lono, Yun & Mihos 2004). In 
the first case, the gas forms shocked, dense filaments which dissi¬ 
pate energy and flow towards the centre where they form a nuclear 
ring which is stabilised by the gravitational potential of the bulge 
- slowing the rate of star formation. In the second scenario, the 
merger between two gas-rich but bulge-less galaxies, the gas forms 
an extended (~ 6 kpc), gravitationally-significant bar-like structure 
which allows the gas to funnel towards the centre. However, with¬ 
out the stabilising influence of the bulges in the progenitor galaxies, 
the gas cannot form a stable ring or disk and its density is further 
increased through the formation of bar-like structures, leading to a 
wide-spread and vigorous starburst. 

These two theoretical scenarios suggest a structural difference 
between the progenitor galaxies may be at the heart of the differ¬ 
ences in behaviour between ULIRGs and SMGs. The suggested 
structural differences are consistent with the expectations of the 
typical galaxies involved in mergers at low- and high-redshifts, 
with the former involving disk-galaxies with significant bulge com¬ 
ponents (e.g. Lilly et al. 1998), while the progenitor galaxies 
of SMGs at high redshift are much more likely to be gas-rich, 
disk-dominated systems with little or no bulge component (Wyse, 
Gilmore & Franx 1997; Ravindranath et al. 2004). If such struc¬ 
tural differences between local and distant starbursts are real, the 
next question is: at what redshift does the transition occur? In il4.3l 


we found that SMGs at 2 ~ 3 have a higher CO detection rate 
than SMGs at 2 ~ 2 - consistent with the former being more gas- 
rich. Secondly, the slight preponderance, albeit tentative, of double- 
peaked CO line profiles at 2 ~ 2 is consistent with that of local 
ULIRGs and in contrast to SMGs at 2 ~ 3 (see II4.3l and II4.4> . Al¬ 
though, these trends are tentative and need a larger survey to con¬ 
firm or disprove them, together they may indicate a difference in 
the physical processes responsible for triggering intense star for¬ 
mation in massive galaxies at 2 ~ 2 compared to 2 ~ 3, with 
major mergers being responsible for much of this activity at lower 
redshifts (as is the case for local ULIRGs), but a separate mode 
(requiring less intense perturbations) capable of triggering similar 
bursts in the more gas-rich systems present at even higher redshifts. 

In I I5.3l we showed that luminous SMGs extend the Leo — Lp,R 
relation of local ULIRGs to higher redshifts and luminosities, and 
accordingly have higher star-formation efficiencies than ULIRGs. 
Although, the star-formation efficiencies found for the SMGs could 
be severely overestimated if an AGN contributes significantly to the 
far-infrared luminosity, the detection of large amounts of molecular 
gas in SMGs along with recent X-ray (Alexander et al. 2003, 2005) 
and radio studies (Chapman et al. 2004) strongly suggest that the 
bulk of the far-infrared emission from SMGs is powered by a large- 
scale starburst and not from an AGN. While CO is a good indicator 
of the total metal-rich H 2 gas reservoir, it may be a worse indi¬ 
cator of the amount of dense gas present (n > 10® cm“®), that 
actually fuels star formation (Carilli et al. 2004). The latter could 
be particularly true in the tidally disrupted giant molecular clouds 
(GMCs) expected in ULIRGs, where a diffuse phase may dominate 
the CO emission but has little to do with star formation (Downes & 
Solomon 1998; Sakamoto et al. 1999). Such a diffuse phase could 
be even more pronounced in SMGs with their more extended dis¬ 
tributions (Chapman et al. 2003c; Small et al. 2004). This would 
explain why the Lfir/Lcq ratio is found to be such a strong function 
of Lfir, increasing for merging systems usually associated with the 
highest far-infrared luminosities. Interestingly, recent work shows 
that the SFE of dense gas, parametrised by the Lfir/Lhcn(i.o) ratio 
(the HCN J = 1 — 0 critical density is 2 x 10® cm“®), remains con¬ 
stant from GMCs all the way to ULIRG system (Gao & Solomon 
2003; Solomon et al. 2003; Carilli et al. 2004). 

The median star-formation rate of the 18 SMGs is (SFR) ~ 
700 Mq yr“^, where we have assumed a Salpeter IMF and a very 
conservative limit corresponding to the starburst contributing only 
50 per cent of the far-infrared luminosity (see Omont et al. 2001, 
but c.f. Alexander et al. 2005). Our findings in 85. H and il5.2l have 
shown that SMGs are massive galaxies with enough molecular gas 
in them to sustain such a large star-formation rate for tsmg ~ 
M{U 2 )/SFR ~ 3 X 10^° Mq /700Mq yr"i ~ 40Myr. By the 
end of such a burst most of the stellar mass corresponding to that 
of a massive spheroid would be in place. A gas depletion timescale 
of tsmg ~ 40Myr is comparable to the typical starburst ages 
(~ 100-200 Myr) derived from photometric modeling of the broad 
band optical/near-infrared colours of SMGs (Small et al. 2004), al¬ 
though we note that large uncertainties are associated with both 
methods. For example, the CO observations only probe the gas 
within the central 10 kpc, and as a result the neutral gas - which 
is likely to be distributed on larger scales - has not been included 
in our estimate of the gas consumption time scale. If the HI gas is 
brought in from > 10 kpc radius at ~ 200kms“^ it would reach 
the central regions in less than 50 Myrs, where it could help sustain 
the vigorous starformation. Furthermore, the above gas consump- 
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tion time scale assumes a continuous starburst until all the gas is 
used - an unlikely scenario since it ignores the negative feedback 
effects from new-born massive stars and supernovae. For example, 
one could imagine the starburst terminating prematurely if the gas 
is removed by starburst- and/or AGN-driven winds only to fall back 
onto the galaxy at a later stage to fuel a second starburst, thus mak¬ 
ing the true gas exhaustion time scale longer. The 40Myr should 
therefore be considered a strict lower limit on the starburst phase of 
SMGs. 

An alternative estimate of the duration of the SMG phase can 
be made from the recent findings by Page et al. (2004) that the 
15 per cent of QSOs at 2 ~ 2 which show absorption in their 
X-ray spectrum are detectable in the submm. Since only 3 per 
cent of radio-identified SMGs are QSOs (Chapman et al. 2005) 
this means that SMGs would have a typical lifetime which is 
~ 0.15/0.03 = 5 times longer than that of QSOs. Adopting a 
QSO lifetime of 40 Myr (Martini & Weinberg 2001) yields a SMG 
life expectancy of tsmg ~ 200 Myr — again larger than the gas 
consumption time scale. 

If we assume that a large fraction of the gas mass of SMGs is 
eventually converted into stars (perhaps through repeated cycles of 
expulsion, infall and star formation), what are the resulting stellar 
masses of the descendents? Clearly SMGs already contain (per¬ 
haps substantial) stellar populations (Smail et al. 2004), however 
the stellar masses for SMGs are difficult to estimate due to sig¬ 
nificant dust extinction and the resulting uncertainties from the de¬ 
generacies between dust reddening and age in evolutionary spectral 
synthesis models. The best estimates to date, based on UK pho¬ 
tometry of 96 SMG and optically-faint /rjy radio galaxies (Smail et 
al. 2004), suggest typical stellar masses of Mstars — 3 x 10^° Mq . 
Thus, combining the gas and stellar mass estimates we find that on 
average SMGs have baryonic masses of ^ 6 x 10^'’ Mq , which is 
comparable to the masses of local early type L* galaxies. 

The discussion above suggests that SMGs are in fact the pro¬ 
genitors of massive spheroids in the present-day Universe and that 
the build-up of the stellar population occurs rapidly (40 g tsmg ^ 
100 Myr), consistent with the old homogeneous stellar populations 
of local ellipticals. Further evidence in support of this is a) the 
strong clustering claimed for SMGs (Blain et al. 2004b), and b) the 
fact that the rest-frame optical properties of SMGs match the bright 
end of the luminosity function of spheroidal galaxies in nearby 
clusters (Smail et al. 2004). 


6 IMPLICATIONS FOR STRUCTURE-FORMATION 
MODELS 

Semi-analytical models of galaxy formation and evolution (Cole 
et al. 1994; Kauffmann et al. 1999; Cole et al. 2000; Somerville 
et al. 2001) have enjoyed a fair degree of success in reproducing 
the properties of galaxies in the local Universe, e.g. the luminosity 
function, the distribution of colours and disk scale lengths of galax¬ 
ies, the observed mix of morphologies and the Tully-Fisher relation 
(e.g. Kauffmann et al. 1999; Cole et al. 2000). However, the key ob¬ 
servable with which to benchmark models of galaxy formation and 
evolution is the mass assembly of galaxies, and in particular the 
assembly of baryonic mass, as a function of redshift. 

Genzel et al. (2003) pointed out that if the extremely large 
baryonic mass (S; 10^^ Mq ) of SMM J02399—0136 was represen¬ 
tative of the bright SMG population, then the high surface density 
of such sources would imply that the abundance of very massive 
baryonic systems at high redshift was about an order of magnitude 
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Figure 9. The co-moving number density of galaxies with baryonic masses 
~ 6 X 10^^ Mq as derived from CO observations of SMGs. The GAL- 
FORM model (Cole et al. 2001; Baugh et al. 2004) predictions of the abun¬ 
dances of galaxies with baryonic masses > 5x, 7x, and 10 X 10^® Mq 
as a function of redshift are shown as thin, medium and thick solid lines, 
respectively. The dotted curve represents the total baryonic matter content 
available in > 10^^ Mq dark matter halos, and is obtained by scaling the 
abundance of halos with the cosmological baryon-to-dark matter density 
(Qi,/Qom = 0.13). This provides a strict upper limit on the number den¬ 
sity of massive baryonic galaxies at a given redshift. The observations can 
be reconciled within the CDM framework, provided that ~ 10 per cent of 
all baryons within dark matter halos are rapidly assembled into galaxies. 


larger than predicted by semi-analytical models (see also Tecza et 
al. 2004). 

Here we repeat this analysis using the 10 SMGs in our sam¬ 
ple which lie in the redshift range 2 = 2-3.5 and have reliable 
gas mass estimates, to put lower limits on the co-moving number 
density of massive galaxies at high redshifts. The redshift interval 
considered translates into ~ 1.5 Gyr in terms of elapsed cosmic 
time, and corresponds to a co-moving volume of 1.8 x 10^ Mpc® 
per square degree. The average submm flux of the 10 SMGs is 
(Sgsojim) = 8 ± 4mJy. The surface density on the sky of SMGs 
with fluxes > 8mJy is 2851^49 deg~^ (Borys et al. 2003). Tak¬ 
ing into account that a) about 60 per cent of the bright SMG pop¬ 
ulation lie within 2 = 2-3.5 (Chapman et al. 2003b, Chapman 
et al. 2005), and b) the CO detection fraction of SMGs in the 
redshift range 2 < 2 < 3.5 is 10/14 (~ 71 per cent), we es¬ 
timate that the co-moving number density of galaxies with bary¬ 
onic masses ^ 6 x 10^° Mq in the redshift interval 2 = 2-3.5 is 
~ S.ljtg ® X 10“® Mpc“®. The errors are estimated by propagating 
the 1-CT limits on the submm flux and number counts through the 
same calculation. 

Since we only observe SMGs during their submm-luminous 
phase, we have to correct the derived space density by a factor 
corresponding to the ratio between the 1.5 Gyr which has elapsed 
over the redshift range 2 = 2.0-3.5 and the typical duration of the 
’SMG phase’. The latter is uncertain, but as we saw in 115.41 it is 
likely to lie in the range 40-200 Myr. Adopting 200 Myr as as a 
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conservative upper value for the submm luminous phase, we esti¬ 
mate a correction factor of ~ 8 . 

In Fig. 1^ we have plotted our estimate of the co-moving num¬ 
ber density of ^ 6 x 10^° Mq systems at 2 ~ 2.8 and the (con¬ 
servative) correction factor which must be applied. Our estimate of 
the abundances of massive galaxies at 2 ~ 2.8 is in good agree¬ 
ment with independent measurements at similar redshifts (Daddi et 
al. 2004; Saracco et al. 2004), and confirms the slow decline in the 
space density of massive galaxies as a function of redshift (Genzel 
et al. 2003; Glazebrook et al. 2004). 

The abundance tracks of > 5 x 10^° and > 7 x 10^'^ Mq sys¬ 
tems as predicted by the most recent GALFORM models (Cole et 
al. 2001; Baugh et al. 2004) are seen to envelope our estimated 
volume density. There are two reasons for the closer agreement 
between the observations and models than previous work: a) the 
typical baryonic mass of bright SMGs is not 2-3 m* but rather 
~ 0.6 m* and b) the new GALFORM model, which employs a top- 
heavy initial mass function (IMF) in order to account for the 850- 
^m number counts (Baugh et al. 2004), also predicts higher abun¬ 
dances of massive galaxies at high redshifts than previous models. 
Other semi-analytical models (Fontana et al. 2004; Granato et al. 
2004) have also modified their recipes so that their models fit the 
data. 

We stress, however, that by using a SMG timescale of 
200 Myr, and not 40 Myr as suggested by the CO observations, our 
estimate of the correction factor is a conservative one. Furthermore, 
if the population of OFRGs, whose faint levels of submm flux is be¬ 
lieved to be due their higher dust temperatures (Blain et al. 2004a; 
Chapman et al. 2005), have similar gas masses as classical SMGs 
(Swinbank et al. 2004; Small et al. 2004), this would double our 
current estimate of the abundance of massive baryonic galaxies at 
high redshifts. 

We conclude that given the small number of sources and crude 
redshift bins used in our calculations, and the large uncertainty on 
the correction factors which have to be applied, there is no evidence 
for a severe discrepancy between observations and the latest theo¬ 
retical predictions. An increase in our sample size would not only 
allow us to confirm this but also allow us to split the sample into 
mass bins and compare with different mass track predictions. 


7 CONCLUSIONS 

We present results from a PdBI CO Survey of SMGs with radio 
counterparts. We successfully detect CO in 4 out of 10 SMGs 
observed, which brings the total number of SMGs detected by 
our survey to 7 out of 13 observed. In addition, we have also 
presented a detection of CO based on archival observations of 
SMM J02396—0134. Combining these 14 sources with four CO- 
detected SMGs from the literature we have compiled a sample of 
18 SMGs observed in CO of which 12 are detected. We use this 
unique sample to derive the bulk gas properties and masses of the 
most luminous SMGs. 

• We find that the SMGs in our sample have a median CO lu¬ 
minosity of (Leo) = (3-8 i 2.0) X 10^° Kkms“^ pc^. This cor¬ 
responds to a molecular gas mass of (M(H 2 )) = (3.0 ± 1.6) x 
10^^ Mq (within R < 2kpc), assuming a conversion factor of 
Xco= 0.8 (Kkm s“^ pc^)“^M 0 and optically thick, thermalised 
line ratios. Although considerable uncertainty is associated with 
Xco it is clear that bright SMGs are amongst the most gas-rich sys¬ 
tems in the Universe. Comparing with local ULIRGs, we find that 
SMGs have molecular gas reservoirs on average about four times 


greater than even the most CO-luminous ULIRGs. We argue that 
this is largely due to the evolution in the molecular gas content of 
the most far-infrared luminous galaxies with redshift. 

• In general, the SMGs in our sample have extremely broad line 
profiles. The median FWHM is 780 ± 320 kms“^ , which is ^ 3 
times larger than the average CO(l—0) line width of local ULIRGs. 
The large line widths, and in some cases multiple-peaked CO spec¬ 
tra, together with the vast amounts of molecular gas, suggests that 
the brightest SMGs are merger events. We argue that the observed 
gas properties are difficult to reconcile with the scenario witnessed 
in local ULIRGs where the bulk of the molecular gas resides in a 
compact circumnuclear disk. 

• We find that the median dynamical mass of the SMG sample 
is {Mdyn} = (1.2 ± 1.5) X 10^^ Mq within the central ~4kpc. 
This suggests that the brightest SMGs are amongst the most mas¬ 
sive galaxies in the distant Universe, comparable in mass to the 
most extreme HzRGs and QSOs, rather than being high-z repli¬ 
cas of local ULIRGs. Taking into account the stellar mass compo¬ 
nent we estimate that the total baryonic mass content of SMGs is 
^ 6 X 10^° Mq . Thus we conclude that not only are SMGs very 
massive baryonic systems, but the baryons can account for a sub¬ 
stantial fraction of the total mass in the central regions. 

• We have shown that the SMGs exhibit a non-linear correla¬ 
tion between far-infrared and CO luminosity, similar to that ob¬ 
served for local ULIRGs. This not only extends the Lco~Lfir re¬ 
lation to higher luminosities (Lhr ~ 10^^”^"^ Lq ) but also shows 
that it holds at the highest redshifts. The main implication of this 
is that SMGs have higher Lfir/Lcq ratios than ULIRGs, and there¬ 
fore possibly higher star-formation efficiencies. However, a clearer 
picture of the star-formation efficiency of SMGs has to await a fu¬ 
ture systematic survey for high-redshift HCN, which will trace the 
dense star-forming gas. 

• From the inferred molecular gas masses we estimate a typical 
gas consumption timescale of > 40 Myr. Given the large uncer- 
tanties involved, an SMG phase of ~ 40 Myr is roughly consistent 
with the best age estimates of the starbursts in SMGs (Tecza et al. 
2004; Small et al. 2003, 2004), which in turn agree favourably with 
independent estimates of the duration of the SMG phase obtained 
by requiring a scenario in which SMGs last 100 — 200 Myrs be¬ 
fore going through a ~ 10 Myr QSO phase at 2 ~ 2, becoming 
massive evolved galaxies at 2 ~ 1, followed by a period of pas¬ 
sive evolution which sees them ending up as old ^ L* ellipticals 
in the present day. Furthermore, the slightly short gas consumption 
time scales can easily be prolonged by invoking negative feedback 
processes, and may be indicative of such processes playing an im¬ 
portant role in the build-up of massive galaxies. 

• From our observations, we place a lower limit on the co¬ 

moving number density of massive baryonic systems in the red¬ 
shift range 2 = 2-3.5 of n{Mhar > 6 x 10^° Mq ) > 

5.1 X 10“® Mpc“® in agreement with results from recent near- 
infrared/spectroscopic surveys (e.g. Saracco et al. 2004; Glaze- 
brook et al. 2004). Given the substantial uncertainties involved, we 
find no significant discrepancy between the data and the predicted 
abundances of massive galaxies at high redshifts. 
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